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Abstract 

Much  of  the  rough  surface  scattering  theory  developed  to  date  considers  only  the 
effects  of  fully  coherent  and  fully  incoherent  illumination  in  the  formation  of  solutions — 
a  problem  studied  in  earnest  since  the  late  1800’s.  In  response,  this  dissertation  extends 
the  theory  currently  available  in  modeling  rough  surface  scattering  to  include  the  effects 
of  partially  coherent  illumination.  Such  illumination  plays  a  pivotal  role  in  our 
understanding  of  active-illumination  systems,  similar  to  those  found  in  directed-energy 
and  remote-sensing  applications,  which  use  the  light  scattered  from  distant  targets  for 
tactical  purposes.  Specifically,  this  dissertation  uses  the  physical  optics  approximation 
(Kirchhoff  boundary  conditions)  to  determine  a  3D  vector  solution  for  the  far-field 
scattering  of  electromagnetic  beam  illumination  with  partial  spatial  coherence  from 
statistically  rough  surfaces.  The  analysis  considers  three  different  material  substrates: 
dielectrics,  conductors,  and  a  perfect  electrical  conductor.  It  also  makes  use  of  a 
Gaussian  Schell-model  form  for  the  incident-field  cross-spectral  density  matrix.  In  so 
doing,  this  dissertation  develops  closed-form  expressions  for  the  scattered  field  cross- 
spectral  density  matrix  with  two  analytical  forms — one  applicable  to  smooth-to- 
moderately  rough  surfaces  and  the  other  applicable  to  very  rough  surfaces.  The  analysis 
shows  that  these  closed-form  expressions  are,  in  general,  complicated  functions  of  both 
the  source  (size  and  coherence  properties)  and  surface  parameters  (surface  height 
standard  deviation  and  correlation  length).  Under  appropriate  conditions,  the  analysis 
also  compares  the  3D  vector  solution  to  previously  validated  solutions  and  empirical 
measurements.  The  results  show  good  agreement. 
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THE  SCATTERING  OE  PARTIALEY  COHERENT  EEECTROMAGNETIC  BEAM 


lEEUMINATION  EROM  STATISTICALEY  ROUGH  SURE  ACES 

1  Introduction 

In  1960,  Theodore  Maiman  invented  the  first  working  laser  [1],  which  originally 
boasted  only  a  few  milliwatts  of  power.  Yet,  by  the  1970s,  laser  powers  reached  the 
megawatt  level  and  the  directed-energy  (DE)  research  community  came  to  life  [2],  The 
technology  found  in  DE  applications  is  inspiring  in  that  it  presents  game-changing 
capabilities  by  offering  systems  with  varying  lethality,  speed-of-light  delivery,  and 
unparalleled  precision  [3-5].  With  this  in  mind,  the  analysis  presented  in  this  dissertation 
hopes  to  aid  the  burgeoning  DE  research  community  and  bring  us  one  step  closer  to 
fielding  an  operational  system  [6,  7]. 

1.1  Problem  statement 

When  using  active-illumination  systems,  more  often  than  not  a  highly  coherent 
laser  beam  propagates  from  the  source  through  the  atmosphere  resulting  in  partially 
coherent  beam  illumination  on  the  target.  This  topic  plays  a  key  role  in  DE  and  remote¬ 
sensing  applications  which  use  the  light  scattered  from  distant  targets  for  tactical 
purposes  [8,9].  Interestingly  enough,  not  much  literature  exists  pertaining  to  the 
scattering  of  partially  coherent  light  from  rough  surfaces. 

In  an  effort  to  bridge  this  gap,  recent  publications  derived  a  2D  scalar-equivalent 
solution  for  the  scattering  of  partially  coherent  beams  from  statistically  rough  surfaces 
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using  the  physical-optics  (PO)  approximation  (Kirchhoff  boundary  conditions)  [10,  11]. 
Specifically,  the  analysis  made  use  of  a  Gaussian  Schell-model  (GSM)  form  in  creating 
the  incident  field  cross-spectral  density  function  (CSDF).  This  allows  one  to  vary  the 
size  and  spatial  coherence  properties  of  the  incident  radiation.  In  so  doing,  the  analysis 
developed  closed-form  expressions  for  the  scattered  field  CSDF  to  observe  the  size  and 
spatial  coherence  properties  of  the  scattered  radiation  in  the  far  field.  The  analysis  also 
validated  these  analytical  expressions  through  computational  simulations  and  showed 
good  agreement  between  the  theoretical  predictions  and  the  numerical  results. 

While  the  2D  scalar-equivalent  solution  is  a  convenient  tool  for  gaining  insight 
into  rough  surface  scattering,  a  complete  understanding  of  the  problem  requires  a  3D 
vector  solution.  With  that  said,  this  dissertation  makes  use  of  the  PO  approximation  to 
determine  a  3D  vector  solution  for  the  far-field  scattering  of  electromagnetic  beam 
illumination  with  partial  spatial  coherence  from  statistically  rough  surfaces.  By 
formulating  the  analysis  in  a  manner  consistent  with  Wolfs  unified  theory  of  coherence 
and  polarization  [12,  13],  all  physical  implications  inherent  in  Wolfs  work  apply  here. 

The  3D  vector  solution  developed  in  this  dissertation  considers  three  different 
material  substrates;  dielectrics,  conductors,  and  a  perfect  electrical  conductor.  In 
addition,  it  uses  a  GSM  form  in  creating  the  incident  field  cross-spectral  density  matrix 
(CSDM).  This  allows  for  the  formulation  of  closed-form  expressions  for  the  scattered 
field  CSDM.  As  such,  the  analysis  shows  that  two  analytical  forms  result  for  the 
scattered  field  CDSM — one  applicable  to  smooth-to-moderately  rough  surfaces  and  the 
other  applicable  to  very  rough  surfaces. 
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Both  analytical  forms  for  the  scattered  field  CSDM  eontain  eomplicated  functions 
of  the  source  parameters  (size  and  eoherence  properties)  and  the  surfaee  parameters 
(surface  height  standard  deviation  and  eorrelation  length).  In  particular,  the  elosed-form 
expression  applicable  to  smooth-to-moderately  rough  surfaces  is  expressed  as  an  infinite 
series.  This  infinite  series  is  slowly  convergent;  however,  one  ean  still  gather  physical 
intuition  from  its  analytical  form.  On  the  other  hand,  the  closed-form  expression 
applicable  to  very  rough  surfaces  is  incredibly  physical,  and  under  certain  cireumstanees, 
maintains  a  GSM  form.  Based  on  these  cireumstanees,  the  analysis  develops  closed-form 
expressions  for  the  angular  speetral  degree  of  coherence  (SDoC)  and  speetral  density 
(SD)  radii.  These  analytical  expressions  also  contain  complicated  functions  of  both  the 
souree  and  the  surface  parameters.  The  analysis  demonstrates  that  for  many  scenarios  of 
interest,  one  ean  approximate/simplify  the  SDoC  radius  as  a  function  of  just  the  source 
parameters  and  the  SD  radius  as  a  funetion  of  just  the  surface  parameters. 

1.2  Dissertation  overview 

Chapters  2  and  3  of  this  dissertation  provide  baekground  information  in  the  form 
of  theory  and  literature  reviews,  respectively.  The  goal  here  is  to  provide  future  researeh 
efforts  with  a  thorough  investigation  of  the  problem  at  hand.  Chapter  4  provides  the 
methodology  used  to  obtain  the  3D  vector  solution  proposed  above.  Here,  the  analysis 
states  all  simplifying  assumptions  and  explains  their  physical  implications.  Chapter  5 
provides  an  exploration  of  the  3D  vector  solution.  The  analysis  given  here  visually 
demonstrates  aspeets  of  the  elosed-form  expressions  and  shows  that  the  results  are 
eonsistent  with  previously  validated  solutions  and  empirical  measurements.  Chapter  6 
provides  a  conclusion  for  this  dissertation  with  a  roadmap  for  future  researeh  efforts. 
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2  Background  electromagnetic  and  optics  theory  review 

Electromagnetic  theory  and  subsequently  optics  theory  involves  the  application  of 
Maxwell’s  equations  to  the  geometry  of  a  specific  problem.  As  such,  the  solutions 
obtained  from  Maxwell’s  equations  determine  the  electromagnetic  vector  fields  present 
within  a  region  of  interest.  With  these  vector  fields,  one  can  then  calculate  quantities  of 
importance  to  their  work.  The  following  chapter  reviews  the  necessary  electromagnetic 
and  optics  theory  needed  to  undertake  the  problem  proposed  above  in  Chapter  1 . 

2.1  Surface  equivalence 

There  are  many  approaches  to  solving  Maxwell’s  equations  for  the  vector  fields 
present  in  an  electromagnetic  scattering  problem.  One  robust  approach  uses  surface 
equivalence.  According  to  Balanis  [14],  Schelkunoff  was  the  first  to  introduce  surface 
equivalence  in  1936  [15].  In  essence,  surface  equivalence  is  a  more  rigorous  extension  of 
Huygen’s  principle  [16],  which,  according  to  Hecht  [17],  states  that  “Every  point  on  a 
propagating  wavefront  serves  as  the  source  of  spherical  secondary  wavelets,  such  that  the 
wavefront  at  some  later  time  is  the  envelope  of  these  wavelets.”  With  this  in  mind, 
surface  equivalence  effectively  defines  equivalent  sources  in  the  form  of  surface  current 
densities  on  a  scattering  object.  Appropriately  defined,  these  surface  current  densities,  in 
addition  to  other  sources,  replicate  the  fields  present  in  a  region  of  interest. 

To  make  this  concept  manifest,  first  consider  the  generic  electromagnetic 
scattering  problem  depicted  in  Eigure  1  [18].  As  shown,  a  primary  source  with  current 
densities,  3’’"  and  ,  radiates  incident  fields,  E'"'’  and  H'"%  which  propagate  in  free 
space  with  index  of  refraction  %=  \  and  impedance  7]^  =  ,  where  £q  and  are 
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the  free-space  permittivity  and  permeability,  respeetively'.  These  known  incident  fields 
illuminate  a  homogeneous,  penetrable  scattering  object  of  volume  V  with  index  of 
refraction  n  and  impedance  rj  .  This  illumination  produces  secondary  sources  in  the 

form  of  current  densities,  and  ,  which  radiate  scattered  fields,  E*'"'  and  H*"' .  It 
is  important  to  note  that  these  scattered  fields  are  unknown.  Furthermore,  a  closed 
surface  S  with  an  outward  pointing  unit-normal  vector  ii  creates  both  an  interior  and 
exterior  region  of  interest.  The  interior  region  contains  the  fields,  E“'  and  H'”' ,  found 
inside  the  closed  surface  S ,  whereas  the  exterior  region  contains  the  fields,  E“^  and 

,  found  outside  the  closed  surface  S  .  These  two  regions  of  interest  create  an  interior 
and  exterior  problem  when  using  surface  equivalence. 


1 .  The  analysis  presented  throughout  this  dissertation  uses  the  MKS  system  of  units  so 
that  fg  =  8.854x10 (farads  per  meter)  and  //g  =  4.;rxl0~^  (henries  per  meter)  [14].  In 
addition,  the  analysis  uses  the  engineering  sign  convention  for  the  time-harmonic 
variations,  i.e.,  u(r,t)  =  U(r,ry)  exp  (yrttt)  =  U(r,ry)e^'** ,  where  U(r,ry)  is  a  position  r 

and  angular  frequency  co  dependent  vector  field  of  interest  within  the  analysis.  Note  that 
sometimes  the  analysis  omits  the  r  and  co  dependence  in  writing  the  vector  fields.  This 
is  done  for  brevity  in  the  notation.  Also  note  that  formulations  which  use  the  physics 
sign  convention,  i.e.,  u(r,^)  =  U(r,ry)exp(-zryt)  =  U  (r,ry)e ,  relate  to  this  work  by  a 

complex  conjugate,  where  j  =  i  =  ^/^ . 

2.  Within  a  homogeneous  space,  the  index  of  refraction  n  and  impedance  rj  relate  to  the 

permittivity  e  and  permeability  ji ,  where  n  =  and  rj  =  [14]. 
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E“'  ^0 
H“'  ^0 


Figure  1.  A  generic  electromagnetic  scattering  problem.  Here,  a  closed  surface 
surrounds  a  scattering  object  creating  an  interior  and  exterior  region  of  interest. 

2.1.1  Interior  problem 

Figure  2  describes  the  interior  problem  when  using  surface  equivalence  [18].  As 
shown,  electric  and  magnetic  surface  current  densities,  and  M"' ,  exist  just  on  the 
inside  of  the  closed  surface  S  and  satisfy  the  following  definitions  [14,  19]; 


4”'  =-hxH'"' 

(1) 

and 

=hxE''''f 

(2) 

These  equivalent  sources  radiate  in  the  absence  of  the  primary  source  and  the  exterior 
fields,  i.e.,  =  0 .  Since  the  exterior  fields  equate  to  zero,  this  allows  the  interior 

region  to  extend  throughout  the  exterior  region  creating  an  unbounded  homogeneous 
space  with  index  of  refraction  n  and  impedance  rj .  In  addition,  the  interior  fields,  E'"' 
and  ,  satisfy  Maxwell’s  equations,  such  that  [19] 
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and 


(4) 


H'"'  =  V  X  A'"^  -j^Ukl+VV]  F'"' , 
nk^  ^ 

where  A!"‘  and  F“'  are  the  interior  magnetie  and  eleetrie  veetor  potentials,  respeetively, 
kf^  =  Inj \  is  the  free-space  wavenumber,  and  is  the  free-space  wavelength^.  Since 

J“'  and  reside  on  a  closed  surface  S  ,  the  vector  potentials  in  Eqs.  (3)  and  (4) 
satisfy  the  following  convolution  integrals  [19]; 


A“=A“(r)  =  j|pr(r'>3(r;r')*' 

(5) 

and 

^in,  ^  ^in,  ^  (r')G(r;r')  , 

(6) 

where  r'  is  the  source  vector,  r  is  the  observation  vector,  and  G(r;r')  is  the  unbounded 
Green’s  function,  such  that 


exp(-7n^o|r-r'|) 

4.;r|r-r  | 

As  a  result,  the  interior  problem  involves  equivalent  sources,  and  ,  which 
replicate  the  interior  fields,  E'"'  and  H'"' ,  in  an  unbounded  homogeneous  space. 


(V) 


3.  Free-space  wavelengths  in  the  optical  regime  typically  range  from  Aq=10^^  pm  in  the 
extreme  ultraviolet  to  Aq=  300  pm  in  the  far  infrared  [20];  thus,  the  analysis  presented 
throughout  this  dissertation  assumes  that  the  free-space  wavenumber  is  much,  much 
greater  than  one,  k^  =  iTrj » 1 . 
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Figure  2.  The  interior  problem  associated  with  using  surface  equivalence.  Here,  the 
exterior  region  contains  null  fields;  consequently,  equivalent  sources  in  the  form  of 
surface  current  densities  replicate  the  interior  fields.  These  equivalent  sources 
effectively  radiate  in  an  infinite  homogeneous  space. 

2.1.2  Exterior  problem 

Similar  to  the  analysis  presented  for  the  interior  problem.  Figure  3  deseribes  the 
exterior  problem  using  surfaee  equivalenee  [18].  Here,  the  eleetrie  and  magnetie  eurrent 
densities,  and  ,  exist  just  on  the  outside  of  the  elosed  surfaee  and  satisfy  the 
following  relationships  [14,  19]; 


Jf  =hxff"' 

(8) 

and 

Mf  =-nxE“'. 

(9) 

These  equivalent  sourees  radiate  in  the  presenee  of  the  primary  souree  while  the  fields  in 
the  interior  region  are  nulled,  i.e.,  E'"'  =  H'"'  =  0  .  Sinee  the  interior  fields  equate  to  zero, 
this  allows  the  exterior  region  to  extend  throughout  the  interior  region  ereating  an 
unbounded  free  spaee.  Moreover,  the  exterior  fields,  and  ,  satisfy  the  following 
superposition  relationships: 


(10) 

and 

jjexf  _  jjmc  _|_  jjs'cf 

(11) 

8 


Solving  Maxwell’s  equations  for  the  unknown  scattered  fields,  and  H*'’' ,  in  Eqs.  (10) 
and  (11),  provides  [19] 


§. 

II 

S' 

1 

<1 

X 

1 

o  |o 

o  to 

+ 

<1 

<1 

(12) 

and 

k  ^ 

to"+VV■)F“^ 

(13) 

0 


where  A“'  and  F“'  are  the  exterior  magnetic  and  electric  vector  potentials,  respectively. 
Since  and  M“'  reside  on  a  closed  surface  S ,  the  vector  potentials  in  Eqs.  (12)  and 
(13)  satisfy  the  following  convolution  integrals  [19]; 


A“=A“(r)  =  j£jr(r')G„(r;r')*' 

(14) 

and 

^  ^  jj^Mf  (r')Go  {Y-y)ds  , 

(15) 

where  here,  (r;r')  is  the  free-space  Green’s  function,  such  that 


exp(-A|r-r1) 


Consequently,  the  exterior  problem  involves  equivalent  sources,  and  ,  which 


radiate  the  scattered  fields,  and  ,  in  the  presence  of  the  known  incident  fields, 
E'"”  and  .  Together  these  fields  replicate  the  exterior  fields,  E*’’"'  and  ,  which 
propagate  throughout  free  space. 
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Figure  3.  The  exterior  problem  associated  with  using  surface  equivalence.  Here, 
the  interior  region  contains  null  fields;  thus,  equivalent  sources  in  the  form  of 
surface  current  densities  radiate  the  scattered  fields.  These  equivalent  sources 
radiate  in  the  presence  of  the  primary  source  and  incident  fields.  Superposition  of 
the  incident  and  scattered  fields  replicate  the  exterior  fields  which  propagate 
throughout  free  space. 

2.1.3  Specialization  to  a  perfect  electrical  conductor 

Figure  4  describes  the  use  of  surface  equivalence  with  a  perfect  electrical 
conductor  (PEC)  [18].  Within  a  perfectly  conducting  material  the  interior  fields  vanish, 
,  SO  that  only  an  exterior  problem  exists  [19].  In  general,  a  PEC  has 
infinite  conductivity,  (7  =  °°  ,  and  the  tangential  electric  field  goes  to  zero  all  along  its 
surface,  i.e.,  hxE“^  =  0  [14].  Based  on  this  knowledge,  the  relationships  given  in  Eqs. 
(8)  and  (9)  simplify  to  the  following  expressions  for  a  PEC: 


Jf  =hxH“' 

(17) 

and 

II 

o 

(18) 

Eurthermore,  Eqs.  (12)  and  (13)  simplify,  such  that 


=  E“"  -  7  ^  (ko'  +  VV  ■ ) 

K 

(19) 

and 

=H''""+VxA“C 

(20) 
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Provided  Eqs.  (17)-(20),  only  an  equivalent  source  Jf‘  is  needed  to  radiate  the  scattered 

fields,  E*""'  and  .  Superposition  with  the  known  incident  fields,  E'"^'  and  H'”'" , 
replicates  the  exterior  fields,  E“'  and  ,  which  propagate  throughout  free  space. 


Figure  4.  The  interior  region  within  a  perfect  electrical  conductor  contains  null 
fields.  As  such,  only  an  exterior  prohlem  exists  when  using  surface  equivalence. 

The  resulting  electric  current  density  radiates  the  scattered  fields  and  superposition 
with  the  known  incident  fields  replicates  the  exterior  fields  which  propagate 
throughout  free  space. 

2.2  Integral  equations 

The  continuity  of  the  tangential  fields  at  the  interface  between  the  exterior  and 
interior  regions  dictates  that  [19] 

h  X  (e“'  -  E'"' )  =  0  ^  h  X  E“'  =  h  X  E'"'  ^  Jf  =  - =r;  (21) 

and  nx(H“'-H“')  =  0^nxH“'=nxH'"'^Mf  =-Mf  =M;^  (22) 

As  such,  in  an  electromagnetic  scattering  problem  using  surface  equivalence,  and 
M^'*'  readily  become  the  primary  unknowns  and  integral  equations  result.  Numerical 
techniques  help  in  solving  these  integral  equations  for  the  unknown  equivalent  surface 
current  densities,  and  M^'*' . 
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2.2.1  Electric-field  integral  equations 

Provided  the  eontinuity  relationship  found  in  Eq.  (21),  the  tangential  eomponents 
of  Eqs.  (3)  and  (12)  simplify,  sueh  that 


nx 


- V  X  -  7  ^  ( yto  +  V V  ■ )  A" 


=  hxE'' 


(23) 


and 


M^^-hx 


-  V  X  -  7  ^  ( n'yto  +  V  V  ■ )  A'"' 


nk. 


=  0 


(24) 


These  eoupled  integro-differential  equations  serve  as  the  eleetrie-field  integral  equations 
(EEIEs)  for  the  unknown  equivalent  surfaee  eurrent  densities,  J^'*'  and  Mj' .  Together, 

Eqs.  (23)  and  (24)  represent  a  linear  system  of  equations  with  two  equations  and  two 
unknowns.  It  is  important  to  note  that  for  a  PEC,  the  tangential  eleetrie  field  goes  to  zero 
all  along  its  surfaee,  i.e.,  nxE“'  =  0 ;  thus,  the  following  EEIE  results  from  Eq.  (19): 


me  ^ 

nxE  =-nx 


(25) 


This  EEIE  is  an  integro-differential  equation  for  the  unknown  equivalent  surfaee  eurrent 


density 


2.2.2  Magnetic-field  integral  equations 

Magnetie-field  integral  equations  (MEIEs)  result  from  the  eontinuity  relationship 
given  in  Eq.  (22).  Speeifieally,  the  tangential  eomponents  of  Eqs.  (4)  and  (13)  simplify, 
so  that 


j;^-nx 


VxA“'-7^(yl:o"+VV-)F‘^ 


nxH" 


(26) 
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and 


(27) 


-J^^-nx  VxA'"^-7^(n\WV-)F'"'  =0. 

L  ^ 

These  coupled  integro-differential  equations  are  the  MFIEs  for  the  unknown  equivalent 
surfaee  eurrent  densities,  and  .  One  should  note  that  Eqs.  (26)  and  (27)  represent 

a  system  of  linear  equations  with  two  equations  and  two  unknowns.  Moreover,  the 
tangential  eomponents  of  Eq.  (20)  ereate  the  following  MFIE  for  a  PEC: 

j;^-nxVxA“'=nxH“C  (28) 

This  MEIE  is  an  integro-differential  equation  for  the  unknown  equivalent  surfaee  eurrent 
density 

2.2.3  Method  of  moments 

The  method  of  moments  (MoM)  is  a  robust  numerieal  approaeh  that  solves  EEIEs 
or  MEIEs  for  the  unknown  equivalent  surface  current  densities,  and  [14,  19,21], 
In  using  the  MoM,  a  series  of  finite  terms  or  basis  funetions  with  unknown  amplitude 
coefficients  effeetively  replace  and  .  This  ereates  a  number  of  algebraic 

expressions  whieh  matrix  algebra  teehniques  readily  solve.  As  a  result,  the  MoM  has  the 
potential  to  formulate  high-fidelity  numerieal  solutions  for  the  unknown  equivalent 
surfaee  eurrent  densities,  3“^  and  . 

2,3  Physical  optics  approximation 

It  is  important  to  remember  that  when  using  surface  equivalenee  in  an 
eleetromagnetie  seattering  problem,  the  equivalent  surfaee  eurrent  densities,  J^'*'  and 

M^'*' ,  radiate  the  unknown  scattered  fields,  E*'"'  and  H*"' ,  whieh  propagate  in  free  space 
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with  index  of  refraction  and  impedance  rj^ .  This  facilitates  the  use  of  a  free-space 


Green’s  function  Gq  (r;r')  [cf  Eq.  (16)];  however,  it  is  still  difficult  to  formulate 
analytical  solutions  for  and  since  they  are,  by  definition,  dependent  on  and 

H*" ,  which  are  unknown.  The  physical-optics  (PO)  approximation  helps  to  alleviate 
these  constraints  [14,  22],  In  essence,  the  PO  approximation  makes  use  of  the 
geometrical-optics  (GO)  approximation  to  formulate  the  current  densities  involved  in  an 
electromagnetic  scattering  problem.  Figure  5  helps  to  further  explain  this  point  [18].  As 
shown,  the  current  densities  formulated  with  the  PO  approximation  equate  to  zero  in  the 
shadow  regions  of  a  scattering  object — an  assumption  which  is  analogous  to  using 
Kirchhoff  boundary  conditions  in  physical  or  wave  optics  [23,  24]. 

When  using  the  PO  approximation  with  surface  equivalence,  one  replaces  the 
unknown  scattered  fields  with  reflected  fields.  Specifically,  ~  and  H*'’'  ~  , 

so  that  the  equivalent  surface  current  densities,  and  become 

=  hxH“^  =  hx(H“"  +  H""' )  -  hx(H''""  +  )  (29) 

and  =  -n  X  E“'  =  -n  x  (e'"'’  +  E""' )  -  -h  x  (e“"  +  ) .  (30) 

These  approximations  assume  that  the  scattering  object  and  its  associated  curvature  are 
large  compared  to  the  wavelength  of  the  incident  fields,  E'”'’  and  Such 

approximations  are  exact  if  the  scattering  object  is  homogeneous,  infinite,  and  planar  [14, 
22].  With  that  said,  the  incident  illumination  follows  the  law  of  reflection — a  direct 

4.  This  is  consistent  with  geometrical  or  ray  optics,  which  emerges  as  the  limit  of 
physical  or  wave  optics  when  the  wavelength  approaches  zero,  ^^0  [20,  22]. 


14 


result  from  the  GO  approximation  [22] — so  that  the  tangential  reflected  fields, 

and  ,  relate  to  the  tangential  incident  fields,  hxE'""  and  hxH'"" ,  in  unique  ways 

for  different  material  substrates,  i.e.,  dielectrics,  conductors,  and  a  PEC.  Appendix  A 


explores  these  relationships  in  more  detail. 


Point 

Source 


Figure  5.  A  macro-scale  description  of  the  physical  optics  (PO)  approximation. 
Here,  a  point  source  illuminates  a  scattering  object  so  that  no  current  densities  exist 
in  the  shadow  region  predicted  hy  the  geometrical  optics  approximation. 


2.3.1  General  problem  for  dielectrics 


As  shown  in  Appendix  A,  the  following  relationships  hold  true  for  dielectrics 


according  to  the  GO  approximation  [18]; 


and 


hxE''‘'^=r5nxE''“  (31) 

hxH''‘^-^=-r5hxH“%  (32) 


where  is  the  known  Fresnel  reflection  coefficient  at  the  closed  surface  5*  of  a  scattering 


object.  Accordingly,  the  equivalent  surface  current  densities,  and  as  given  in 
Eqs.  (29)  and  (30),  simplify  so  that 


(33) 
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(34) 


and  M;" --(l  +  r5)nxE“^ 

This  dictates  that  when  using  the  PO  approximation  for  a  dielectric  material,  the  known 
tangential  incident  fields,  hxE'”'"  and  hxH'”'",  in  addition  to  the  known  Fresnel 
reflection  coefficient,  ,  are  all  that  is  needed  in  determining  the  analytieal  forms  of  J^'*' 

and  M5" . 

2.3.2  Specialization  to  conductors 

For  very  good  conductors,  the  conductivity  approaches  infinity,  cr  ^  00  [14],  As 
such,  the  tangential  incident  electric  field  approximates  to  zero  all  along  the  surface  of  a 
eonducting  material,  i.e.,  nxE'"""  ~  0,  and  Eqs.  (33)  and  (34)  simplify  to  the  following 
relationships: 

J?=(l-f,)nxH'"  (35) 

and  MJ=0.  (36) 

Thus,  the  analysis  simplifies  from  that  of  dielectrics.  Only  the  electric  equivalent  surface 
eurrent  density  3“^  radiates  when  using  the  PO  approximation  for  a  eonducting  material. 


2.3.3  Specialization  to  a  perfect  electrical  conductor 
As  shown  in  Appendix  A,  the  following  relationships  hold  true  for  a  PEC 
aecording  to  the  GO  approximation  [18]: 


nxE''"^  =-nxE“" 

(37) 

and 

nxH^‘'^=nxH“C 

(38) 

Consequently,  the  equivalent  surface  eurrent  densities,  J^'*'  and  ,  as  given  in  Eqs. 
(29)  and  (30),  simply  so  that 


16 


(39) 


j;*'  -  2nxH“" 

and  =  0 .  (40) 

This  dictates  that  when  using  the  PO  approximation  for  a  perfectly  conducting  material, 
the  known  tangential  incident  magnetic  field  hxH'"'"  is  all  that  is  needed  in  determining 
the  analytical  form  of  J^'*' . 

2.4  Mathematical  techniques 

Mathematical  techniques  exist  which  further  simplify  the  analysis  beyond  the  PO 
approximation.  To  help  make  these  mathematical  techniques  unambiguous,  first  consider 
the  3D  electromagnetic  scattering  setup  described  in  Figure  6.  The  analysis  also  refers  to 
this  setup  as  the  macro-scale  scattering  geometry.  Here,  the  vector,  p  =  xx-l-  mu  ,  points 
from  the  source  plane  origin  to  a  transverse  location  in  the  source  s  ,  since  v  =  0  in  the 
source  plane;  the  vector,  f  =  Ox  -  yy  -l-  zz  ,  points  from  the  source  plane  origin  to  the 
surface  plane  origin;  the  vector,  r'  =  x'\  +  y'y  +  z'z  ,  points  from  the  surface  plane  origin 
to  a  point  on  the  closed  surface  S  ;  and  the  vector,  r  =  xx  +  yy  +  zz  ,  points  from  the 
surface  plane  origin  to  an  observation  point.  This  setup  plays  a  pivotal  role  in  employing 
the  far-field  approximation,  plane-wave  spectrum  representation,  and  method  of 
stationary  phase,  all  of  which  are  mathematical  techniques  which  greatly  simplify  the 
analysis  when  used  under  the  right  assumptions. 
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Figure  6.  A  generic  setup  for  3D  electromagnetic  scattering  problems.  Here,  the  x 
axes  align  in  both  the  source  plane  and  the  surface  plane.  This  assumes  isotropy  in 
the  material  substrate. 

2.4.1  Far-field  approximation 

In  the  far  field,  r-r'  is  approximately  parallel  to  r  ,  as  shown  in  Figure  7  [18]. 
Additionally,  |r  -  r'| »  1 ,  and  the  following  approximations  result  [14]: 

I  ,1  - 7  [r  -  f  ■  r'  for  phase  variations 

r-r' =Vr'+/'-2r-r'-^  ^  ,  (41) 

[  r  for  amplitude  variations 


where  r  =  |r| ,  r  =  |r1 ,  and 


f  ■r''  =  x''sin(0)cos(^)  +  ysin(^)sin(^)  +  z''cos(^) . 


(42) 


Provided  Eqs.  (41)  and  (42),  the  free-spaee  Green’s  fimetion  (r;r') ,  as  given  in  Eq. 
(16),  simplifies  so  that 


C.  (r;  O  =  • 

4.7rr-r  ''  — 


Ajtr 


(43) 
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From  Eq.  (43),  the  unknown  scattered  fields,  E*""'  and  ,  as  formulated  above  in  Eqs. 
(12)  and  (13),  satisfy  the  following  approximations  [14]; 


=-VxF“'-7^(yi:o +VV-)A“' 

-wry  set  ,  TP  SCt 

(44) 

and 

=  V X  -j^[kl+VV) = 

K 

1  TT-SC? 

(45) 

where 

(46) 

H-=^fxE;;^ 

V, 

(47) 

Hr=-/o^xp(-7V)...^^ 

t/o  Anr  ^ 

(48) 

and 

Ef  =-;7ofxHf . 

(49) 

These  expressions  depend  on  the  far-field-exterior  magnetic  and  electric  vector 
potentials,  N“'  and  ,  such  that 


N“  =  N“  (r)  =  117?  (r')exp(;V  ■  r')*' 

(50) 

and 

L“  =  L”  ( r )  =  £  M?  (r')  exp  {j\r  r')ds' . 

(51) 

It  is  important  to  remember  that  the  PO  approximation  simplifies  the  analysis  so  that  the 
known  incident  fields,  E'"'"  and  H'"'' ,  help  in  determining  the  analytical  form  of  the 
equivalent  surface  current  densities,  and  in  Eqs.  (50)  and  (51).  Moving 

forward,  one  can  use  the  plane -wave  spectrum  representation  to  account  for  E'""  and 


z 


Figure  7.  The  micro-scale  geometry  describing  far-field  electromagnetic  scattering. 

Here,  the  observation  point  extends  well  past  what  is  depicted  and  only  a  small 
portion  of  the  closed  surface  is  displayed.  This  corresponds  to  a  zoomed-in 
description  of  Figure  6. 

2.4.2  Plane-wave  spectrum  representation 

In  order  to  determine  the  analytieal  form  of  the  equivalent  surfaee  eurrent 
densities,  and  M^'*' ,  the  analysis  must  first  aeeount  for  the  known  ineident  fields,  E'"'" 


and  H'”'" ,  whieh  propagate  from  the  souree  s  and  illuminate  the  elosed  surfaee  S  [ef. 
Figure  6].  With  this  in  mind,  one  ean  write  the  ineident  eleetrie  field  E'"'"  in  terms  of  its 
speetrum  using  the  plane  wave  speetrum  representation  [25].  The  following 
expressions  result  using  the  maero-seale  seattering  geometry  deseribed  in  Figure  6; 


^inc  ^  ^inc  ^  f  ?  r^inc  \  '  (f  +  F  )]  dk‘^^ V  >  0  (52) 


and  T:“  = 


e  e 


[kf , c )  =  n  (p)  [i'  +  Cu )]  dxdu,  V  >  0  .  (53) 
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Here,  k'“  =  ytf  x  +  ytf  u  +  is  the  ineident  propagation  veetor  and  v  >  0  is  the  souree- 

free  half  spaee  with  index  of  refraetion  Hq  =  1  and  impedanee  .  Sinee  the 

divergenee  of  the  ineident  eleetrie  field  equals  zero  in  a  souree-free  half  spaee,  i.e., 

V  ■  E'""^  =  0 ,  it  follows  that  in  the  spatial-frequency  domain. 


rjimc  j^inc rji inc  _|_  j^incrj-^inc  _|_  j^incrpinc  q 


(54) 


This  physically  states  that  the  spectrum  of  the  incident  electric  field  Vf"  is  perpendicular 


to  the  incident  propagation  vector  k'"'" .  Thus,  for  ^  0 , 


7  inc  1  inc 

rj-^inc  _  rj-iinc  rj^inc 


(55) 


so  that  the  x  and  u  components  of  T"‘'  uniquely  provide  the  v  component. 

In  a  similar  fashion,  the  following  relationships  provide  the  plane-wave  spectrum 
representation  for  the  incident  magnetic  field  H'"'"  and  its  spectrum  [25]; 


H 


^H“(r)  =  ^i  iTr(r.C)Mp[-lk'”'(f  +  r)]rf<:r<^C'.  ''SO  (56) 


and  \  \n‘’'^{p)Qx^\^j[k^x  +  kl"%)\dxdu,  v>0.  (57) 


From  Maxwell’s  equations,  VxE'”'’  =  -yryo/ZoH'"'' ,  so  that  in  the  spatial-frequency 
domain, 

k'”'’xTf  =  ryo//oTf .  (58) 

Here,  rq,  =  2.;rVQ  is  the  free-space  angular  frequency  and  Vq  is  the  free-space  frequency. 

In  an  equivalent  form,  k'"'’  =  k^k'"''  =  CDQ^J~JU^k‘'“' ,  where  k'"'’  is  the  incident  unit- 
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propagation  vector.  Consequently,  the  spectrum  of  the  incident  magnetic  field  T^' 
relates  to  the  spectrum  of  the  incident  electric  field  in  the  following  ways: 


^  rjiZMC 


(0,11^ 


**-k*”^  X 


^  rjimc 


% 


(59) 


It  follows  that  is  also  perpendicular  to  the  incident  propagation  vector  k'"'" ;  namely, 


j^mc  rpmc 


1  incrpinc  .  i  incrpinc  .  i  incrpinc  _ ^ 


(60) 


SO  that  for  kT  ^  0 . 


"/iv 


(61) 


This  says  that  both  the  x  and  u  components  of  uniquely  provide  the  v  component. 

Before  moving  on  in  the  analysis,  it  is  important  to  note  that  the  expressions  given 
in  Eqs.  (53)  and  (57)  for  the  incident  spectrums,  and  ,  are  mathematically 
equivalent  to  taking  the  two-dimensional  Fourier  transform  of  the  incident  fields  in  the 
source  plane,  i.e.,  E'"'"  (p)  and  H'""^  (p) ,  and  going  to  the  spatial-frequency  domain. 

With  that  said,  the  expressions  given  in  Eqs.  (52)  and  (57)  allow  one  to  then  determine 
incident  fields  at  any  observation  point,  i.e.,  E'""  (r)  and  H'""  (r) .  Physically,  this  is 

analogous  to  summing  up  the  contributions  of  a  bunch  of  forward  propagating  plane 
waves  which  originate  from  the  source  plane  [20].  It  is  also  important  to  note  that  the 
plane- wave  spectrum  representation  often  results  in  rather  complex  integral  expressions. 
In  practice,  one  must  employ  additional  mathematical  techniques,  such  as  the  method  of 
stationary  phase,  to  solve  these  complex  integral  expressions. 
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2. 4. 3  Method  of  stationary  phase 

In  using  the  plane-wave  spectrum  representation  (along  with  the  far-field  and  PO 
approximations),  integrals  of  the  following  form  often  result  [18]; 

b 

P(k)  =  \f(x)  exp  [  jkg  (x)]  dx,  as  A:  ^  oo  .  (62) 


Here,  / (x)  is  slowly  varying  in  the  interval  \a,b\  and  Ag-(x)  is  rapidly  oscillating 
except  near  special  points  where  the  rate  of  change  of  g(x)  is  stationary  within  the 
interval,  i.e.,  where 


^g{x)  =  g{x)  =  i).  (63) 

dx 

These  special  points  are  called  critical  points  of  the  first  kind  [26].  Away  from  these 
points,  kg  (x)  is  rapidly  oscillating  and  the  positive  and  negative  contributions  of  the 

integrand  effectively  cancel  out.  In  this  case,  an  asymptotic  mathematical  technique 
known  as  the  method  of  stationary  phase  helps  in  solving  the  integral  formulated  in  Eq. 
(62). 


The  initial  analysis  assumes  that  there  is  only  one  critical  point  of  the  first  kind; 
namely,  at  x  =  Xg  and  that  /  (x)  and  g  (x)  are  both  continuous  and  well  behaved  in  the 

interval  [a,h] .  Subsequently,  the  following  conditions  must  hold  true:  g(xQ)  0, 
g'(xo )  =  0 ,  and  g" [x^^]  0  ,  so  that  upon  expanding  in  a  Taylor  series. 


g(x)  =  g(xo)  +  g  (Xo)(x-Xo)  +  g''(Xo)^^  +...^g(Xo)  +  g''(Xo)^^  (64) 
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and  /(x)  =  /(xo)  +  /'(xo)(x-Xo)  +  /''(xo)^^^^  +  ...-/(xo),  (65) 

since  /  (x)  is  slowly  varying.  Substituting  Eqs.  (64)  and  (65)  into  Eq.  (62),  the  method 
of  stationary  phase  dictates  that  [26] 

F{k)^F^^\k)  +  F^^\k),  as^^oo.  (66) 

Here,  [k)  =  f{x^) exp \_jkg (xq )]  j  exp  jk  g" {x^ )  dx  (67) 


is  the  contribution  from  the  eritieal  point  of  the  first  kind  at  x  =  Xg , 


and  (k)  is  the 


contribution  from  the  end  points.  These  end  points  are  called  critical  points  of  the  second 
kind  [26];  however,  the  present  analysis  negleets  to  formulate  their  eontributions^.  With 
some  mathematical  prowess,  Eq.  (66)  evaluates  to  the  following  expression  [26]; 


F{k)^F^'\k)  = 


(68) 


where 


ifg''(xo)>0 

ifg''(xo)<0’ 


(69) 


In  general,  if  there  are  multiple  critical  points  of  the  first  kind  present  in  the  analysis,  then 
their  individual  eontributions  sum  together. 

The  analysis  leading  up  to  Eqs.  (68)  and  (69)  assumed  one-dimensional 
integration;  however,  the  method  of  stationary  phase  extends  to  n-dimensional  integration 


5.  Critical  points  of  the  seeond  kind  do  not  eome  into  play  beeause  of  the  nature  of  the 
fields  assumed  in  this  research  effort. 
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[27,  28].  As  such,  the  contribution  [k)  from  the  critical  point  of  the  first  kind  at 


x  =  Xo  =(xio,X2o,...,x„o)  becomes 


^  In 

V  k  , 


njl 


/(^o) 


Det{a'g(xo)} 


1/2 


exp[7lg(i 


,)]exp  7-sgn{a^g(xo)}  ,  (70) 

V  4  / 


9'g(x) 

dxjdxj 

3^g(x) 

where  a2g(xj  = 

dxjdx^ 

(71) 

9'g(x) 

9"g(x) 

dx„dx2 

X=Xo 

Det{---|  denotes  the  determinant  operation,  and  sgnj^l}  = denotes  the 
signature  of  a  real  symmetric  non-degenerate  matrix  A  .  Here,  /l+  [A]  are  the  number  of 

positive  and  negative  eigenvalues  of  A  . 

2,5  Coherence 

The  fields  of  interest  in  electromagnetic  scattering  problems  are  often  random  in 
nature.  Goodman  refers  to  such  fields  as  optical  disturbances  [29].  Of  primary  concern 
in  the  statistical  analysis  of  optical  disturbances  is  coherence.  In  essence,  coherence 
describes  the  degree  to  which  one  point  in  a  given  optical  disturbance  relates  to  any  other 
point  within  the  optical  disturbance  in  time  or  space.  An  optical  disturbance  is  coherent 
when  there  is  a  fixed  relation  between  one  point  and  all  other  points  within  the  optical 
disturbance.  On  the  other  hand,  an  optical  disturbance  is  then  incoherent  when  there  is 
no  fixed  relation  between  one  point  and  any  other  point.  Statistical  properties  that  fall 
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somewhere  between  the  preceding  descriptions  provide  for  a  partially  coherent  optical 
disturbance.  Mathematically,  one  realizes  coherence  through  correlation  functions 
r  (rj ,  Tj ;  )  •  These  correlation  functions  are,  in  general,  dependent  on  two  points  in 

space,  Tj  and  ,  or  two  instances  of  time,  h  and  • 

2.5.1  Self-coherence  function 

When  analyzing  temporal  coherence,  an  individual  uses  what  Goodman  refers  to 
as  the  self-coherence  function  r(r,T)  [29].  Specifically, 

j  r/2 

r(r,r)  =  lim—  [  u{r,t  +  T)u  {r,t)dt  =  {u{r,t  +  T)u  {r,t)\ ,  (72) 

^  -Tjl 

which  is  simply  the  time  autocorrelation  of  an  analytic  function  u  (r,t)  at  a  single  point 

in  space  r .  Throughout  the  analysis,  u  (r,t)  represents  the  optical  disturbance  of  interest 

and  u  (r,t)  represents  the  complex  conjugate  of  that  optical  disturbance^.  Note  that  in 

writing  Eq.  (72),  one  assumes  that  the  optical  disturbance  is  emanating  from  a  point 
source,  so  that  only  temporal  coherence  effects  plays  a  role.  As  such,  the  temporal 
quantity,  r  =  -  h  ,  is  the  time  interval  of  interest  in  quantifying  temporal  coherence.  Put 

simply,  the  self-coherence  function  r(r,T)  gives  a  distinct  gauge  for  temporal  coherence 

provided  r  <K  ~  Itt! So) ,  where  is  the  coherence  time  and  Ary  is  the  finite  angular 

bandwidth  of  the  optical  disturbance  [29].  The  process  used  to  measure  temporal 
coherence  helps  in  explaining  this  point  further. 

6.  Note  that  the  scalar  field  analysis  presented  here  holds  for  vector  fields,  i.e.,  each 
component  of  the  vector  field. 


26 


In  practice,  coherence  measurements  require  the  interference  of  light  using  optical 
devices  called  interferometers.  The  type  of  interferometer  used  depends  highly  on  the 
type  of  coherence  to  be  measured  for  a  given  optical  source.  For  example,  when 
temporal  coherence  is  of  concern,  light  from  a  point  source  is  interfered  with  a  delayed 
version  of  itself  This  type  of  interference  requires  amplitude  splitting  of  the  light.  A 
Michelson  interferometer  achieves  this  type  of  interference  and  is  readily  described 
throughout  the  optics  literature — Goodman’s  treatment  is  particularly  insightful  [29].  In 
the  detection  plane  of  a  Michelson  interferometer,  the  irradiance  /  (r)  scales  with  the 

self-coherence  function  r(r,T) .  This  says  that  one  can  measure  temporal  coherence 

through  the  interferenee  of  light. 

2.5.2  Mutual-coherence  function 

When  analyzing  spatial  coherence,  an  individual  uses  what  Wolf,  Goodman,  and 
many  others  refer  to  as  the  mutual  coherence  function  (MCF)  r(rj,r2,T)  [13,  29].  In 
particular, 

T{r^,r^,T)  =  {u{r^,t  +  T)u  (73) 

which  is  a  time  cross  correlation  of  an  analytic  function  u  (r,t)  at  two  points  in  space,  Tj 
and  Tj .  When  dealing  with  a  single  point  in  space  r  ,  Eq.  (73)  reduces  to  a  self¬ 
coherence  function  F (r, r) ,  as  given  in  Eq.  (72),  so  that  in  general,  the  MCE  E (rj ,  Tj , r) 

is  more  robust  in  quantifying  coherence.  In  writing  Eq.  (73)  and  similarly  Eq.  (72),  one 
assumes  that  the  optical  disturbance  is  statistically  stationary,  at  least  in  the  wide  sense. 
This  means  that  the  average  optieal  disturbance  has  no  explieit  time  dependence;  instead. 
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the  MCF  r(rj,r2,T)  depends  only  on  time  difference,  ^  =  not  the  actual  values  of 
tj  and  ^2  •  Physically,  this  is  analogous  to  steady-state/continuous-wave  operation  of  the 
optical  source  [13]. 

In  writing  Eq.  (73),  one  also  assumes  that  the  optical  disturbance  is  emanating 
from  an  extended  source,  so  that  both  temporal  and  spatial  coherence  effects  play  a  role. 
Temporal  coherence  effects  play  a  role  in  the  definition  of  the  MCF  E  (rj,r2,T)  because 

there  is  the  potential  for  optical  path-length  differences  between  the  extended  source  and 
the  two  points  in  space,  Tj  and  r2 .  These  optical  path-length  differences  are  negligible 

when  there  is  symmetry  between  the  extended  source  and  the  two  points,  Tj  and  r2 ,  and 
when  the  light  is  quasimonochromatic  or  narrowband,  such  that  A<w«c  d) ,  where  (d  is  the 
mean  angular  frequency  of  the  optical  disturbance  [13,  29].  When  these  conditions  are 
met,  the  analysis  treats  the  temporal  properties  within  the  MCF  r(rj,r2,T)  separately, 
viz., 

r  (ri ,  r2 ,  t)  -  J  (ri ,  r2 )  exp  {jar)  (74) 

and  J(rj,r2)  =  r(rj,r2,T  =  0)  =  ^M(rj,t)M*(r2,t)^.  (75) 

As  such,  the  mutual  intensity  J(rj,r2 )  gives  a  distinct  gauge  for  spatial  coherence 
provided  the  two  points  in  space,  Tj  and  r2 ,  situate  themselves  within  the  spatial 
coherence  area,  where  X  is  the  mean  wavelength  of  the  optical 

disturbance  and  is  the  solid  angle  subtended  from  the  extended  source  to  the  two 
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points  [13,  29].  The  process  used  to  measure  spatial  coherence  helps  in  explaining  this 
point  further. 

When  spatial  coherence  is  of  concern,  one  would  want  to  interfere  the  light  from 
an  extended  source  with  a  spatially  shifted,  but  not  delayed  version  of  itself  [29].  This 
type  of  interference  requires  wavefront  splitting  at  two  separate  points.  The  Young’s 
double  slit  experiment  achieves  this  type  of  interference  and  is  readily  described 
throughout  the  optics  literature — the  treatments  of  Goodman  and  Wolf  are  particularly 
insightful  [13,  29].  In  the  detection  plane  of  Young’s  double  slit  experiment,  the 
irradiance  scales  with  the  MCF  r(rj,r2,T) ,  and  if  the  setup  allows  for  it,  the 

irradiance  /(r)  scales  with  the  mutual  intensity  J(rj,r2) .  This  says  that  an  individual 

can  measure  spatial  coherence  through  the  interference  of  light. 

2. 5. 3  Complex  degree  of  coherence 

Normalizing  the  MCF  r(rj,r2,T) ,  as  given  in  Eq.  (73),  an  individual  obtains  a 


quantity  referred  to  as  the  complex  degree  of  coherence  (CDoC)  y[r^,r2,T) ,  where 


r(rpr2,T) 


_ _ 

Vr(ri,ri,T  =  0)r(r2,r2,T  =  0) 


(76) 


Note  that  the  complex  degree  of  (self)  coherence  7(r,r)  follows  from  Eq.  (72)  when 
dealing  with  a  single  point  in  space  r  [29],  and  similarly,  the  (equal-time)  complex 
degree  of  coherence  7  (ri,r2 )  follows  from  Eq.  (75)  when  dealing  with  symmetry  in  the 
optical  setup  and  narrowband  light  [13].  Eurthermore,  one  can  relate  the  visibility 


29 


V  of  the  irradiance  ,  detected  in  their  respective  interferometer,  to  the  CDoC 


7(rj,r2,T)  using  the  following  relationship  [13,  29]; 


Vir  =  r  t)| 

max[/(r,T)]  +  min[/(r,T)] 


(77) 


Together,  Eqs.  (76)  and  (77)  say  that  the  magnitude  of  the  CDoC  |y(rj,r2,T)|  provides  a 

normalized  unit  of  measure  for  the  amount  of  coherence  (temporal  or  spatial)  in  an 
optical  disturbance  at  two  points  in  space,  Tj  and  Tj  ,  and  some  time  difference  T  = 

For  example,  if  |y(rj,r2,T)|  =  1 ,  two  different  points  in  space  are  correlated  and  the 

optical  disturbance  is  fully  coherent;  however,  if  |y(rj,r2,T)|  =  0  ,  two  different  points  in 
space  are  uncorrelated  and  the  optical  disturbance  is  incoherent.  A  partially  coherent 
optical  disturbance  then  satisfies  0  <  |y(rj,r2,T)|  <  1 . 

2.5.4  Cross-spectral  density 

The  cross-spectral  density  (CSD)  W  (rj,r2,ty)  is  an  alternative  way  of  analyzing 


spatial  coherence  [13].  Explicitly, 


r  (rj ,  Tj ,  r)  =  ^  f  IE  (rj ,  Tj ,  ^y)  exp  (y^yr)  t/r 


271  { 


(78) 


and 


IE  (rj ,  r2 ,  yt)  =  J  r  (rj ,  r2 ,  r)  exp  {-jcor)  dr , 


(79) 


such  that  the  MCE  E  (tj  ,  r2 ,  r)  and  the  CSD  IE  (tj  ,  Tj  ,  &>)  form  a  Fourier  transform  pair. 
This  says  that  the  CSD  W is  a  way  to  analyze  spatial  coherence  in  the  space- 
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frequency  domain  as  opposed  to  the  space-time  domain  with  the  MCF  r(rj,r2,T) . 
Moreover,  Wolf  derives  the  following  result  [13]; 

W {r^,r^,a))  =  {r^,0})U\r^,0})) ,  (80) 

which  says  that  the  CSD  W  (r^rj,  ttt)  is  the  cross  correlation  function  of  an  ensemble 
\U  (r,ty)|  of  sample  functions  U  {r,0)) .  These  sample  functions  are  the  space-  and 
angular- frequency-dependent  part  of  a  monochromatic  optical  disturbance,  i.e., 
u{Y,t)  =  U  {y,(o)  QX^[jCOt). 

It  is  import  to  remember  that,  in  general,  the  Fourier  transform  of  an  optical 
disturbance  does  not  exist  because  it  is  not  absolutely  integrable.  However,  the  Wiener- 
Khintchine  theorem  states  that  for  a  random  process  that  is  zero  mean  and  at  least  wide- 
sense  stationary,  the  autocorrelation  and  the  spectral  density  form  a  Fourier  transform 
pair  [13,  29].  This  is  an  important  point  in  the  analysis  because  when  dealing  with  a 
single  point  in  space  r ,  Eq.  (80)  reduces  to  an  expression  for  the  spectral  density 
,  where 

S (r, (o)  =  W (r, r,  co)  =  (u (r,  co)U*  (r, co)'^ .  (81) 

Thus,  the  spectral  density  S  (r,  co)  is  a  way  to  analyze  self  coherence  in  the  space- 
frequency  domain  as  opposed  to  the  space-time  domain  with  the  self-coherence  function 
r (r,T) .  The  process  used  to  measure  spatial  coherence  in  the  space-frequency  domain 
helps  in  explaining  this  point  further. 

When  measuring  spatial  coherence  in  the  space-frequency  domain,  one  again  uses 
Young’s  double  slit  experiment.  Narrow-band  filters  placed  behind  the  slits  ensure  that 
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the  optical  disturbances  emanating  are  space-  and  angular-frequency-dependent 
ensembles,  |l7 (rj,&>)|  and  \U {r^,co)^  .  This  allows  an  individual  to  consider  the 
spectrum  of  the  light  in  the  detection  plane  instead  of  irradiance.  Specifically,  one  can 
measure  the  spectral  density  S'  (r,  ry) .  Wolf  shows  that  in  the  detection  plane  of  this 

modified  Young’s  double  slit  experiment  [13],  the  spectral  density  S  {y,co)  scales  with 

the  CSD  W (rj,r2,ty) .  This  says  that  an  individual  can  measure  spatial  coherence 

through  the  interference  of  light  in  the  space-frequency  domain. 

2. 5. 5  Spectral  degree  of  coherence 

Normalizing  the  CSD  W ,  as  given  in  Eq.  (80),  an  individual  obtains  a 
quantity  referred  to  as  the  spectral  degree  of  coherence  (SDoC)  p  ,  where 


W{y„y^,co) 

^W{y^,y^,cu)W{y2,Y2,cu) 


(82) 


As  such,  one  can  then  relate  the  visibility  V  (r,  co)  of  the  spectral  density  S  (r,  co) , 


detected  in  the  modified  Young’s  double  slit  experiment,  to  the  SDoC  ju{y^,Y2,co)  using 


the  following  relationship  [13]; 


max  [5*  (r,  ft))]  -  min  [5*  (r,  ry)] 
max  [iS  (r ,  ry)  J  -l-  min  [5*  (r ,  ry)  J 


jU{Y„Y^,a)) 


(83) 


Provided  Eqs.  (82)  and  (83),  the  magnitude  of  the  SDoC  |^(rj,r2,ty)|  provides  a 

normalized  unit  of  measure  for  the  amount  of  spatial  coherence  in  an  optical  disturbance 
at  two  points  in  space,  Tj  and  r2 ,  and  angular  frequency  w.  Eor  instance,  if 
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|//  (r^Tj,  ^y)|  =  1 ,  two  different  points  in  spaee  are  eorrelated  and  the  optieal  disturbanee  is 

spatially  eoherent;  however,  if  |//(rj,r2,ty)|  =  0 ,  two  different  points  in  spaee  are 
uneorrelated  and  the  optieal  disturbanee  is  spatially  ineoherent.  A  spatially  partially 
eoherent  optieal  disturbanee  then  satisfies  0  <  |//  (Tj  ,  Tj  ,  ty)|  <  1 . 

2. 5. 6  Gaussian  Schell-model  source 

Refereneing  Figure  7,  in  the  souree  plane  at  v  =  0 ,  the  CSD  W {r^,r^,w)  takes  the 
following  form: 

lF(R,p2,^y)  =  (f/(pi,<y)f/*(p2,ftt)),  (84) 

where  Pj2  =Xj2X  +  Mj2U.  Consequently,  the  CSD  W[p^,p2,ci))  of  a  Gaussian  Sehell- 
model  (GSM)  souree  takes  the  following  form  [13]: 


lF(ft,P2,ty)  =  ^5(fi,ry)^5(p2,ty)//(p2-fi,ty), 


(85) 


sueh  that 


S'(p,ty)  =  exp 


(86) 


and 


ju{p,co)  =  Qxp 


(87) 


Note  that  the  parameters  ,  w ,  and  i  are  spaee  independent  but  are,  in  general, 
dependent  on  angular  frequeney  co .  This  dependenee  is  omitted  for  brevity  in  the 
notation.  Also  note  that  upon  substituting  Eq.  (85)  into  Eq.  (82),  the  magnitude  of  the 

SDoC  |^(p2,Pi,ty)|  beeomes 
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(88) 


|//  (Pz ,  ft ,  ^y)|  =  |//  (P2  -  R ,  ^y)|  =  exp 

which  depends  only  on  the  distance  between  two  points  and  not  on  the  points  themselves. 
This  is  the  classic  characteristic  of  a  Schell-model  source  [13]. 

Similar  to  the  Gaussian  laser  beam  source  [20],  the  three  parameters  ,  w ,  and 
i  physically  describe  the  GSM  source.  For  instance,  the  source  beam  width  w  is  the 
radial  distance  |p|  where  the  source  magnitude  falls  to  l/e^  its  initial  on-axis  value. 
This  gives  a  nice  gauge  for  the  physical  size  of  the  emanating  beam.  Likewise,  the 
source  coherence  length  £  is  the  distance  between  two  points  |p2  -  ft  |  where  the 

magnitude  of  the  SDoC  |^(p2,ft,ty)|  falls  to  l/e^  its  initial  on-axis  value.  This  is  a 
direct  result  of  the  relationship  found  in  Eq.  (88).  In  practice,  if  |p2  -  ft  |  <K  f  ,  then  the 
two  points  are  correlated  and  the  GSM  source  is  spatially  coherent;  conversely,  if 
|p2  -  ft  I  »  f  ,  then  the  two  points  are  uncorrelated  and  the  GSM  source  is  spatially 

incoherent.  Partial  spatial  coherence  then  satisfies  0  <  [pz  -  ft  |  <  f . 

In  using  the  GSM  source  formulated  in  Eqs.  (85)-(87),  the  analysis  is  tractable  for 
a  variety  of  fields  of  practical  interest.  Eor  example,  the  GSM  source  reduces  to  a  point 
source  when  the  source  beam  width  approaches  zero,  w  ^  0 ,  or  a  plane  wave  when  the 
source  beam  width  approaches  infinity,  w^o°  .  One  can  also  use  the  GSM  source  to 
model  spatially  coherent  Gaussian  laser  beams.  Here,  an  individual  allows  the  source 
coherence  radius  to  approach  infinity,  £  On  the  other  hand,  when  the  source 

coherence  radius  approaches  zero,  £  ^0 ,  one  obtains  a  spatially  incoherent  Gaussian 


2  P2  Pl 
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beam  source.  This  simplicity  and  versatility  makes  the  GSM  source  ideal  for 
investigations  concerned  with  spatial  coherence. 

2.5. 7  Cross-spectral  density  matrix 

When  analyzing  spatial  coherence  in  the  space-frequency  domain  with 
electromagnetic  vector  fields,  one  uses  the  cross-spectral  density  matrix  (CSDM) 

W  (rj,r2,ty)  [13].  In  general,  the  cross-spectral  density  matrix  CSDM  W(rj,r2,tt>)  is 
the  dyadic  (outer  product)  created  from  electric  field  vectors  of  the  following  form; 


E  (r„  ^y)  =  ,  <y)  X  +  ,  ^y)  u  (/  =  1, 2) 


E,{r^,co) 


(/  =  1.2) 


(89) 


such  that 


W  ( ri ,  r2 ,  yt)  =  (e  ( ri ,  y;)  ( r2 ,  yt) ) 


/r^x(rp«) 


_  {e^{y„co)eI{y„co))  I^E^{y„co)eI{y^,co)) 
_{Eu  (rp  El  (r2,  yt))  (£„  (ri,  y;) El  {y^,co)) 
=  (^E^{y^,(o)eI{y^,(o))  (m  =  x,M;  n  =  x,M) 

=  ^mn(*'p*'2’y^)  (m  =  x,M;  n  =  x,M) 


(90) 


where  t  denotes  Hermitian  conjugate.  In  Eq.  (89),  E^{yj,co)  and  E^{Y^,co)  are  members 

of  statistical  ensembles  which  are  at  least  wide-sense  stationary,  and  in  referencing 
Figure  7,  are  analytic  functions  in  two  mutually  orthogonal  directions  perpendicular  to 
the  direction  of  propagation,  i.e.,  the  v  direction.  This  says  that  the  vector-field  result 
presented  in  Eq.  (90)  is  analogous  to  the  scalar-field  result  given  above  in  Eq.  (80). 
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Accordingly,  the  SD  andtheSDoC  ^(rj,r2,ty)  are  determined  from  the  CSDM 

W (rj,r2,tt>)  using  the  following  relationships  [13]; 


S  (r,  co)  =  Tr  1 W  (r,  r,  ft»)| 

(91) 

and 

Tr{w(r„r„®)} 

p-— - T  - 7’ 

^Tr{w(r„r„®)}^Tr{w  r„r„®)} 

(92) 

where  Tr  {•  •  •)  denotes  the  traee  operation  and  r  =  Tj  2 .  This  says  that  the  veetor-field 
result  presented  in  Eqs.  (91)  and  (92)  directly  relate  to  the  scalar-field  results  given  above 
in  Eqs.  (81)-(83).  The  magnitude  of  the  SDoC  |//(rj,r2,&>)|  resulting  from 
electromagnetic  vector  fields  also  provides  a  normalized  unit  of  measure  for  the  amount 
of  spatial  coherence,  i.e.,  0  <  |//(rj,r2,ty)|  <  1 . 

Referencing  Eigure  7,  the  CSDM  W  (rj,r2,  ttt)  of  a  Gaussian  Sehell-model  (GSM) 
source  takes  the  following  element-based  form  [13]; 

^mn  (pp  p2  .  «)  =  (R  .  (p2  .  (p2  -  Pi  ^  «)  (»!  =  X,  M;  n  =  X,  M  )  ,  (93) 


/^mn(p2-R,«)  =  ^m„eXp 


2p2-R 


(m  =  x,m;  n  =  x,m) 


Note  that  the  element-based  parameters  Al^,  w^,  5^„,and  are  space  independent  but 
are,  in  general,  dependent  on  angular  frequency  a> .  This  dependence  is  omitted  for 
brevity  in  the  notation.  Also  note  that  the  CSDM  W {r^,r2,co)  given  in  Eqs.  (93)-(95)  is 
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analogous  to  the  scalar-field  result  given  above  in  Eqs.  (85)-(87);  however,  there  are 
additional  constraints,  viz., 

5mn  =  1  when  m  =  n  ,  (96) 


S  „  <  1  when  mj^n 


(97) 


and 


B  =B  , 

mn  nm  ’ 


mn  nm 


(98) 


(99) 


Nonetheless,  the  GSM  source  presented  in  Eqs.  (93)-(99)  is  ideal  for  investigations 
concerned  with  spatial  coherence. 

2,6  Polarization 

Given  electromagnetic  vector  fields  and  the  CSDM  W(rj,r2,ty) ,  as  defined 

above  in  Eq.  (90),  polarization  relationships  result.  The  first  polarization  relationship  of 
interest  is  the  space-  and  angular- frequency-dependent  degree  of  polarization  (DoP) 
P[r,a))  [13].  Particularly, 


P{r,Q)) 


i 


^  4Det{w(r,r,^y)} 
(Tr{w(r,r,6;)})' 


(100) 


where  again,  Det{-  •  •)  denotes  the  determinant  operation  and  r  =  Tj  j .  In  general,  the 
DoP  P[r,co)  provides  a  normalized  unit  of  measure  for  the  amount  of  polarization  in  an 
optical  disturbance  [13,  30].  When  P[r,co)  =  \,  the  optical  disturbance  is  polarized, 
whereas  when  P{r,O})  =  0,  the  optical  disturbance  is  unpolarized.  Partial  polarization 
then  satisfies  0<E’(r,ry)<l. 
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The  second  polarization  relationship  of  interest  is  the  space-  and  angular- 
frequency-dependent  angle  of  polarization  (AOP)  y/[Y,co) .  Specifically, 


^^(r,ry)  =  ^tan  ' 


2Re{y„(r,r.a))} 

KAr,r,0))-W,Jr,r,a)) 


,  [-nl2<y/<nl2).  (101) 


This  angle  is  depicted  in  Figure  8  in  terms  of  a  polarization  ellipse.  The  semi-major  and 
semi-minor  axes  of  this  ellipse  satisfy  the  following  relationships  [13]; 


and 


(r,  a)  =  ^1^4  (r ,  r,  (r,  r ,  -  ify Jr,  r,  cd)J 


4 (r, r , ^y)}]  +[w^{r,r,cu)-W^^{r, r, co)\ 

'  (r,  CD)  =  ^|^4|ify„(r,r,^y)|"  +  [r„(r,r,6i)-iF„„(r,r,fti)]' 


4 (r , r, ^y)}]  +[w^{r,r,CD)- (r , r, yi)]' 


(102) 


(103) 


respectively.  As  such,  the  space-  and  angular-frequency-dependent  ellipticity  e(r,CD) 


follows  as 


•  (r,  cd)  =  =  tan  [z{r,  co)~\ ,  (-;r/4  <x<nlA), 

a(r,ft»)  L  V  /j  V 


(104) 


where  j(r,&))  is  the  ellipticity  angle,  which  is  also  depicted  in  Figure  8.  In  general,  the 
ellipticity  e[r,  cd)  provides  a  normalized  unit  of  measure  for  the  polarization  state  of  the 
optical  disturbance  [23,  30].  When  e(r,CD)  =  l,  the  semi-major  and  semi-minor  axes  of 
the  polarization  ellipse  equal  each  other,  i.e.,  a  (r,  CD)  =  b  (r,  cd)  .  This  corresponds  to  a 


circularly  polarized  optical  disturbance.  On  the  other  hand,  when  £•  (r,  ry)  =  0  ,  the  semi- 
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minor  axis  of  the  polarization  ellipse  equals  zero,  i.e.,  b{r,(o)  =  0.  This  eorresponds  to  a 
linearly  polarized  optieal  disturbanee.  An  elliptieally  polarized  optieal  disturbanee  then 
satisfies  0  <  <  1 . 


Figure  8.  A  description  of  the  polarization  geometry  used  in  the  analysis. 

The  third  polarization  relationship  of  interest  is  the  angular- frequeney-dependent 

two-point  Stokes  veetor  s(rj,r2,ty)  [13].  Per  se,  the  eomponents  of  this  veetor  satisfy  the 

following  relationships: 

s,{y„y^,co)  =  W^{y„y^,co)  +  {r„Y^,co) 

s,{Y„Y^,co)  =  W^{r„r^,co)- (r^ , , ty) 
s^{y,,y^,co)  =  W^^{y,,y^,co)  +  W^{y,,y^,co) 
s,{y^,y„cd)  =  j[w^^{y„y„co)- (rj ,  ,  ty)] 

A  (t ,r2, cu)  =  (a,  (ri ,  ty) El  {y^ ,  a))  +  (a„  {y, ,  cd) eI  (r^ ,  a)) 

s,  (ri , ,  ty)  =  (a,  (r^ ,  ty) E]  (r^ ,  ty))  -  (a„  (r^ ,  ty)  El  (r^ ,  ty)) 
or 

^2  (ri ,  ,  ty)  =  {E^  (r^ ,  ty) E]  (r2 ,  ty))  +  (A„  (r^ ,  ty) E]  (r2 ,  ty)) 

^3  (t  ,  r2 ,  ^y)  =  j\{K  (ri  ,co)eI{y^,  co))  -  {e^  {y,  ,  a)  El  {y,  ,  co))_ 
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History  shows  that  the  single  point  Stokes  veetor  s  (r ,  r,  ry) ,  where  r  =  Tj  2 ,  is  a  very 

versatile  tool  in  terms  of  analyzing  polarization  [23,  30].  For  example, 

,  ,  Jsf(r,r,CD)  +  sUr,r,CD)  +  sUr,r,CD) 

P(r,co)=^  1  V  ^  ^  ; — 2  V  ; — ^  ^  \  (107) 

5o(r,r,6t) 

whieh  says  that  one  ean  obtain  the  DoP  P[r,  co)  two  separate  ways  within  the  analysis. 
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3  Background  rough  surface  scattering  literature  review 

As  mentioned  in  Chapter  1 ,  the  purpose  of  this  dissertation  is  to  extend  the  rough 
surfaee  seattering  literature  to  inelude  the  effeets  of  partially  eoherent  eleetromagnetie 
beam  illumination.  In  support,  reeent  publieations  derived  a  2D  sealar-equivalent 
solution  for  the  seattering  of  partially  eoherent  beams  from  statistieally  rough  surfaees 
using  the  physieal-opties  (PO)  approximation  [10,  11].  These  publieations  serve  as  the 
basis  for  this  dissertation;  however,  modern-day  researeh  in  rough  surfaee  seattering 
dates  baek  to  the  work  of  Lord  Rayleigh  around  the  turn  of  the  20**^  eentury  [31-33]. 

With  this  said,  one  ean  distinguish  the  published  literature  in  rough  surfaee  seattering  into 
two  main  eategories.  The  first  eategory  deals  with  the  researeh  predominately  eoneemed 
with  the  seattering  of  fully  eoherent  and  fully  ineoherent  illumination  from  rough 
surfaees,  whereas  the  seeond  eategory  deals  with  the  researeh  predominately  eoneemed 
with  the  seattering  of  partially  eoherent  illumination  from  rough  surfaees. 

3,1  Fully  coherent  and  fully  incoherent  illumination 

Several  different  researeh  eommunities  eome  to  mind  when  reviewing  the  rough 
surfaee  seattering  literature  pertaining  to  fully  eoherent  and  fully  ineoherent  illumination. 
The  first  eouple  identify  themselves  with  the  rough  surfaee  seattering  researeh  performed 
by  the  opties  and  photonies  eommunities  for  metrology  and  manufaeturing  applieations. 
The  text  written  by  Stover  highlights  this  point  [34].  Conversely,  the  seeond  eouple 
identify  themselves  with  the  rough  surfaee  seattering  researeh  performed  by  the  radio- 
frequeney/mierowave  and  visible/near-infrared  eommunities  for  synthetie  aperture  radar 
and  remote  sensing  applieations.  The  three-volume  text  by  Ulaby  et  al.  highlights  this 
point  [35].  With  some  exeeptions,  the  eommon  approaehes  employed  by  these  researeh 
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communities  are  the  geometrieal-opties  (GO),  linear-systems,  perturbation,  PO,  and  full- 
wave  methods.  One  may  refer  to  works  of  Beekmann  and  Spizziehino  [36],  Ishimaru 
[37],  Ogilvy  [38],  Voronovieh  [39],  Wamiek  and  Chew  [40],  Elfouhaily  and  Guerin  [41], 
Nieto- Vesperinas  [42],  Maradudin  [43],  and  Fung  and  Chen  [44],  for  exeellent 
summaries  on  rough  surfaee  seattering  teehniques  using  fully  eoherent  and  fully 
ineoherent  illumination. 

3.1.1  Geometrical-optics  methods 

When  employing  GO  methods,  one  typieally  uses  a  bidireetional  refleetanee 
distribution  funetion  (BRDF)  or  its  polarimetrie  eounterpart,  a  polarimetrie  BRDF 
(pBRDF),  to  model  rough  surfaee  seattering.  Nieodemus  was  the  first  to  introduee  the 
BRDF  in  1965  [45].  Defined  in  radiometrie  terms,  the  BRDF  is  the  refieeted  radianee 
divided  by  the  ineident  irradianee  [34].  As  sueh,  the  BRDF  typieally  eharaeterizes  how 
light  refieets  from  surfaees  in  terms  of  a  speeular  and  diffuse  eomponent  or  a  polarized 
and  unpolarized  eomponent  [46].  Both  empirieal  and  analytieal  BRDFs  exist  in  praetiee. 
Measurements  help  in  formulating  empirieal  BRDFs.  For  example,  in  preparation  for  the 
NASA  Apollo  missions,  the  analysis  of  light  seattered  from  the  lunar  surfaee  led 
researehers  to  eonelude  that  the  moon’s  surfaee  is  eomposed  of  a  partieulate  material 
[47].  On  the  other  hand,  both  PO  and  GO  approximations  help  in  formulating  analytieal 
BRDFs  [22].  Using  the  GO  approximation,  the  seminal  BRDF  paper  is  that  of  Torranee 
and  Sparrow  [48],  whereas  the  oft-refereneed  pBRDF  paper  is  that  of  Priest  and  Meier 
[49].  Many  other  models  exist  based  on  their  work.  These  inelude  BRDFs  for 
applieations  in  passive  visible/near-infrared  remote  sensing  [50]  and  eomputer  graphies 
[51] — Sun’s  literature  review  is  partieularly  insightful  [52]. 
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3.1.2  Linear-systems  methods 

In  the  late  1970s,  Harvey  and  Shack  developed  a  linear-systems  formulation  of 
rough  surface  scattering  based  on  scalar  diffraction  theory  [53-55].  In  this  approach,  a 
surface  transfer  function  characterizes  the  scattering  process  much  like  the  optical 
transfer  function  does  for  aberrations  found  within  an  imaging  system.  The  Fourier 
transform  of  this  surface  transfer  function  then  yields  a  scattered  radiance  distribution 
function  closely  related  to  the  BRDF.  In  the  late  1980’s,  Harvey  et  al.  modified  this 
theory  to  include  the  effects  of  grazing  incidence  at  X-ray  wavelengths  [56].  This  helped 
in  the  design  of  X-ray  telescopes.  Most  recently,  Krywonos  et  al.  modified  the  theory 
once  again  to  a  non-paraxial  regime  [57,  58].  This  scalar  non-paraxial  linear-systems 
formulation  of  rough  surface  scattering  claims  to  produce  accurate  results  for  rougher 
surfaces  than  the  theories  based  on  perturbation  methods  and  for  larger  incident  and 
scattered  angles  than  the  theories  based  on  PO  methods  [59-61]. 

3.1.3  Perturbation  methods 

The  perturbation  approach  to  rough  surface  scattering  models  the  surface 
roughness  as  a  small  perturbation  relative  to  the  case  of  a  perfectly  smooth  surface.  As 
such,  this  approach  requires  that  the  surface  roughness  be  small  compared  to  the 
wavelength  of  the  incident  radiation  [37].  The  literature  credits  Rice  with  the 
groundbreaking  paper  on  this  subject  [62];  however,  it  is  important  to  note  that  Lord 
Rayleigh  initiated  the  use  of  many  of  the  mathematical  techniques  [31-33].  Thus,  the 
literature  often  refers  to  the  perturbation  formulation  of  rough  surface  scattering  as 
Rayleigh-Rice  theory.  It  is  also  important  to  note  that  different  approaches  found  within 
the  literature  tend  to  yield  similar  results  up  to  a  fifth-order  perturbation  expansion  [63]; 
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nonetheless,  perturbation  methods  are  the  oldest  and  most  widely  used  in  the  rough 
surfaee  seattering  literature. 

3.1.4  Physical-optics  methods 

The  PO  approaeh  to  rough  surfaee  seattering  uses  the  PO  approximation  [14,  22], 
whieh  is  analogous  to  using  Kirehhoff  boundary  eonditions  in  physieal  or  wave  opties 
[23,  24],  This  is  done  so  that  instead  of  satisfying  the  exaet  boundary  eonditions,  as  is 
done  with  perturbation  methods,  the  field  and  its  normal  derivative  simplify  on  the 
seattering  surfaee.  Aeeordingly,  this  approaeh  does  not  require  that  the  surfaee  roughness 
be  small  eompared  to  the  wavelength  of  the  ineident  radiation  [35].  The  literature 
typieally  eredits  Beekmann  with  the  trailblazing  work  on  this  subjeet  [36],  and  one  often 
sees  the  title  of  Beekmann-Kirehhoff  theory  or  the  Kirehhoff  approximation  used  in 
praetiee.  It  is  important  to  note  that  the  PO  approximation  typieally  allows  an  individual 
to  ealeulate  elosed-form  expressions  where  other  approximations/theories  would  not. 

Sueh  is  the  ease  when  eonsidering  the  seattering  of  fully  eoherent  laser  beam  illumination 
from  rough  surfaees  [64-68]. 

3. 1. 5  Full-wave  methods 

When  employing  full- wave  methods,  one  typieally  uses  the  method  of  moments 
[14,  19,  21],  the  finite  differenee  time  domain  [19,  69,  70],  or  the  finite  element  method 
[19,  71,  72]  to  satisfy  Maxwell’s  equations  and  model  rough  surfaee  seattering.  This 
problem  has  a  rieh  history  that  dates  baek  to  the  late  1970’s.  Some  of  the  early  notable 
work  in  this  field  is  that  of  Bahar  [73-76],  Axline  and  Fung  [77],  Thorsos  [78],  and  Collin 
[79,  81].  The  topieal  review  written  by  Warniek  and  Chew  outlines  many  sueh  full- wave 
teehniques  [40]. 
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3,2  Partially  coherent  illumination 

The  rough  surface  scattering  literature  pertaining  to  partially  coherent 
illumination  is  the  sole  result  of  the  proliferation  of  laser-based  systems,  such  as  those 
found  in  active-illumination  systems  for  directed-energy  and  remote-sensing  applications 
[82-85].  In  recent  times,  the  statistical  behaviors  of  the  laser-target  interaction;  in 
particular,  the  resulting  speckle  patterns,  gained  considerable  interest.  Since  the  presence 
of  speckle  is  typically  detrimental  in  applications  involving  coherent  light,  techniques  for 
suppressing  speckle  naturally  followed.  Some  of  the  early  notable  literature  in  such 
fields  as  metrology  and  remote  sensing  include  the  research  efforts  of  Dainty  [86],  Fujii 
and  Asakura  [87,  88],  Pedersen  [89],  Goodman  [90],  Parry  [91],  and  Yoshimura  et  al. 
[92].  A  recent  text  written  by  Goodman  reviews  many  such  techniques  [93]. 

One  way  to  suppress  speckle  in  active-illumination  systems  is  to  use  partially 
coherent  light  instead  of  fully  coherent  laser  light.  As  a  result,  this  research  topic  is 
becoming  more  and  more  popular  due  predominately  to  the  work  of  Wolf  in  creating  his 
unified  theory  of  coherence  and  polarization  [12,  13].  This  unified  theory  helps  in 
explaining  correlation-induced  changes  in  coherence,  polarization,  and  spectrum  of 
partially  coherent  light.  In  particular,  much  of  the  published  literature  uses  the  properties 
of  a  partially  coherent  electromagnetic  beam  whose  cross-spectral  density  matrix 
possesses  a  Gaussian  Schell-model  (GSM)  form  [13].  As  the  name  implies,  Schell  was 
the  first  to  conjecture  such  an  electromagnetic  source  in  1961  [94,  95].  Since  then,  much 
effort  has  gone  into  understanding  the  physics  behind  GSM  sources/beams.  Some  of  the 
published  literature  includes  research  in  realizability  conditions  [96-98],  experimental 
generation  [98-100],  numerical  simulation  [102-104],  free-space  propagation  [105-107], 
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turbulent  propagation  [108-1 10],  or  sources/beams  of  similar  form  [111-113].  The 
topieal  review  written  by  Gbur  and  Visser  [114]  and  a  reeent  text  written  by  Korotkova 
[115]  are  most  thorough  in  review. 

In  regards  to  the  seattering  of  partially  eoherent  illumination,  most  of  the  eurrent 
literature  deals  with  the  seattering  from  low-eontrast  surfaees,  i.e.,  where  the  index  of 
refraetion  differs  only  slightly  from  unity  [116-125].  These  are  seattering  surfaees  in 
whieh  the  Born  approximation  is  valid  [13,  23,  1 14,  115] — the  topieal  review  written  by 
Zhao  and  Wang  thoroughly  reviews  this  problem  [126].  In  view  of  this,  there  are  far 
fewer  publieations  on  the  seattering  of  partially  eoherent  illumination  from  rough 
surfaees.  Of  the  published  work  to  date,  the  following  approaehes  are  eommon;  the 
phase-sereen  model,  the  ABCD  matrix  formulation,  and  the  eoherent-mode 
representation. 

3.2.1  Phase-screen  methods 

In  1975,  Goodman  developed  a  phase-sereen  formulation  of  rough  surfaee 
seattering  based  on  sealar  diffraetion  theory  [90,  93].  In  this  approaeh,  a  phase-sereen 
transmittanee  funetion  eharaeterizes  the  seattering  proeess  mueh  like  an  aperture 
transmittanee  funetion  does  in  physieal  or  wave  opties  [24].  Hoover  and  Gamiz  most 
reeently  employed  this  approaeh  [127].  In  so  doing.  Hoover  and  Gamiz  assumed 
idealized  quasimonoehromatie  plane-wave  illumination.  This  allowed  for  the  applieation 
of  the  generalized  Van  Cittert-Zernike  (VCZ)  theorem  to  the  mutual  intensity  funetion  on 
the  phase-sereen  surfaee.  The  VCZ  theorem  relates  the  irradianee  to  the  mutual  intensity 
through  a  Fourier  transform  [13,  29].  Hoover  and  Gamiz ’s  work  ultimately  lead  to  the 
formulation  of  a  generalized  BRDF  solution  whieh  was  the  sum  of  a  eoherent  and 
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incoherent  component;  however,  their  work  did  not  direetly  aecount  for  partially  eoherent 
beam  illumination. 

3.2.2  ABCD-matrix  methods 

The  text  written  by  Andrews  and  Phillips  best  describes  the  principles  behind  the 
ABCD-matrix  approaeh  to  rough  surfaee  seattering  [128].  In  general,  the  ABCD-matrix 
approach  describes  paraxial  wave  propagation  through  any  complex  optical  system. 

When  modeling  the  rough  surfaee  seattering  using  the  ABCD-matrix  approaeh,  a  phase- 
screen  transmittance  function  again  characterizes  the  scattering  process.  However,  the 
inclusion  of  a  soft-Gaussian  aperture  in  the  model  aecounts  for  the  size  of  the  scattering 
surface  and  accompanying  diffraction  effects.  Korotkova  discusses  this  point  in  her  texts 
[115,  129].  The  researeh  of  Hansen  et  al.  [130]  and  Yura  and  Hanson  [131]  used  this 
approach  to  look  at  rough  surface  scattering  from  a  target  which  produced  partially 
developed  speekle.  Wu  and  Cai  also  deseribed  an  approaeh  to  sensing  the  seatter  from 
rough  surfaces  using  ABCD-matrix  methods  and  partially  eoherent  beam  illumination  via 
the  GSM  formulation  [132];  however,  this  work  is  only  applicable  to  small-angle 
seattering  geometries  with  very  rough  surfaces  [133-135]. 

3.2.3  Coherent-mode  methods 

The  text  written  by  Ostrovsky  best  deseribes  the  prineiples  behind  the  eoherent- 
mode  approach  to  rough  surface  scattering  [136].  Huttunen  et  al.  used  this  approach 
along  with  the  PO  approximation  to  look  at  the  seattering  from  two-dimensional  miero- 
struetured  media  [137],  i.e.,  an  isolated  groove  or  slit  in  a  perfectly  conducting  material 
substrate.  This  unique  approach  to  rough  surface  scattering  may  prove  useful  for  future 
research  efforts. 
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4  Methodology  for  the  3D  vector  solution 

Figure  1  describes  the  geometry  used  to  obtain  a  3D  vector  solution  for  the 
problem  proposed  above  in  Chapter  1 .  As  shown,  a  zero  mean  2D  sample  function 
h  =  h[x,y)  describes  the  surface  height  at  the  rough  interface  with  standard  deviation 
and  correlation  length  f  ^ .  This  gives  rise  to  a  statistically  rough  surface  S  .  Spatially 
partially  coherent  electromagnetic  beam  illumination  (parameters  given  below)  emanates 
from  the  source  plane  specified  by  the  coordinates  [x,u,v) ,  which  are  different  from  the 

surface-plane  coordinates  (x,y,z).  As  such,  the  vector,  p  =  xx -I- mu  ,  points  from  the 

source  plane  origin  to  a  transverse  beam  location  since  v  =  0  in  the  source  plane;  the 
vector,  =  Ox-y^y  +  z^z ,  points  from  the  source-plane  origin  to  the  surface-plane 

origin;  and  the  vector,  r  =  xx  +  yy  +  zz  ,  points  from  the  surface  plane  origin  to  an 
observation  point.  Note  that  in  the  source  and  the  surface  planes,  the  x  axes  align,  which 
assumes  that  the  surface  of  the  homogeneous  medium  is  statistically  isotropic  [38,  43] 
with  impedance  rj .  Above,  the  medium  is  free  space  with  impedance  //g . 

4,1  Incident  field  cross-spectral  density  matrix 

As  mentioned  above,  spatially  partially  coherent  electromagnetic  beam 
illumination  emanates  from  the  source  plane.  With  this  in  mind,  the  analysis  uses  a 
Gaussian  Schell-model  (GSM)  form  for  the  incident  field  cross-spectral  density  matrix 

(CSDM)  W'  (pi,p2)  [13,  94],  such  that 
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where  (•  •  •)  denotes  eorrelation,  t  denotes  Hermitian  eonjugate,  and  *  denotes  eomplex 


eonjugate.  In  Eq.  (108),  the  element-based  parameters  A^  and  A^  are  the  beam 


amplitudes  in  the  x  and  u  directions,  respectively,  is  the  source  width,  and  the 
element-based  parameters  and  are  the  correlation  amplitude  and 


correlation  length,  respectively.  Note  that  B^^  follows  additional  constraints  [13,  97]; 
namely, 

^mn  =  1  when  m  =  n 

|5jj^|<l  whenm?^:n.  (109) 


Also  note  that,  in  general,  the  incident  field  E'  and  the  parameters  A^,  w^,  and  B^ 

are  radian  frequency  O)  dependent  [13];  however,  the  analysis  omits  this  dependence  for 
brevity  in  the  notation^. 


7.  The  analysis  presented  in  this  chapter  uses  the  MKS  system  of  units  in  addition  to  the 
engineering  sign  convention  for  the  time-harmonic  variations  (cf.  Footnote  1,  p.  5).  In 
addition,  some  of  the  notation  is  simplified  from  that  presented  in  Chapter  2,  e.g., 
gmc  _  Pqj.  brevity  in  the  notation. 
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(a) 


(b) 

Figure  9.  The  macro-scale  (a)  and  micro-scale  (b)  scattering  geometry  of  a  2D 
statistically  rough  surface  S  of  length  2L  and  width  2L. 


50 


4,2  Scattered  field 


The  present  analysis  uses  the  PO  approximation  to  develop  a  far-field  expression 
for  the  scattered  field  E*  (r) .  For  this  purpose,  one  can  write  the  incident  field  E'  (r)  in 

terms  of  its  spectrum  T'  {k[,k[ )  using  the  plane -wave  spectrum  representation  [25]. 

Using  the  macro-scale  scattering  geometry  given  in  Figure  9a,  the  following  expressions 
result: 

r(r)  =  ^i  i r(k‘,k'y'^'-<^>^'dk‘dK  (110) 

(2;r)  LL 

and  T'  [k[X)=]]  E'  (p)  ,  (111) 


which  are  valid  in  the  source-free  half  space  where  v  >  0  .  In  Eqs.  (110)  and  (111), 
k'  =  k^lC  =  k[x  +  k'Ji  +  A:' V  is  the  incident  propagation  vector,  k^  =  Inj \  is  the  free- 
space  wavenumber,  and  is  the  free-space  wavelength. 

For  most  directed-energy  and  remote-sensing  engagement  scenarios,  all  of  the 
observation  points  of  interest  are  in  the  far  field.  As  such,  the  scattered  electric  field 
E*  (r )  depends  on  the  far-field  vector  potentials,  L  (r)  and  N  (r ) ,  using  the  following 

o 

relationships  : 


(112) 


8.  In  the  far  field,  r  >  2D^  j ,  where  D  ~2L  and  r:^L-,  consequently,  the  analysis 
neglects  all  contributions  to  the  scattered  field  E^  that  are  in  the  radial  f  direction 
because  their  contributions  scale  as  l/ r" ,  where  n  =  2, 3, . . . ,  and  are  negligible  [14]. 
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L(r)  =  J|M(r')e^''»"'''j/, 

5 

(113) 

and 

(114) 

s 


In  Eqs.  (1 12)-(1 14),  0  and  are  unit  vectors  in  the  polar  (vertical  polarization)  and 
azimuth  (horizontal  polarization)  directions,  respectively,  M(r')  and  are  the 
equivalent  surface  current  densities,  respectively,  and  the  vector,  r'  =  x'x  +  yy  +  h'z  , 
points  from  the  surface-plane  origin  to  a  point  on  the  statistically  rough  surface  S  . 

Using  the  micro-scale  scattering  geometry  given  in  Figure  9b  and  the  PO 


approximation  [14,  22], 

M(r')--M'-E‘'(r')  (115) 

and  J(r')-— J'-E‘(r'),  (116) 

Vo 

where  J'  and  M'  are  dyadics,  such  that 

J'  =  h'x  s'(l-?]')p'-p'|l-r/ js'  (117) 

and  M'  =  h'x  s'(l-l-r[)s'-l-p'(l-l-?]|')p'  .  (118) 


In  Eqs.  (117)  and(118),  n  =  n[x',y')  is  the  2D  unit  outward  normal  vector  given  by 


-hA-hyy  +  i 

Jl  +  h^,  +  hy, 


(119) 


where 


dh'  _  dh  {x',  y')  _  dh'  _  dh  {x',  y') 


(120) 
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Furthermore,  s'  and  p'  are  the  unit  perpendieular  and  parallel  veetors,  whereas  and 

rji  are  the  eorresponding  Fresnel  refleetion  eoeffieients,  respeetively.  Refereneing  the 
miero-seale  seattering  geometry  in  Figure  9b,  the  following  relationships  result: 


s  =  p  =  sxk'  p''=-sxk''. 

k'xh 


(121) 


Thus,  in  arriving  at  the  relationships  found  in  Eqs.  (117)  and  (118),  one  must  use  the  GO 
approximation  [22];  speoifieally  the  law  of  refleetion,  sueh  that  -h  ■  k'  =  h  ■  k'' 

Based  on  Eq.  (119),  it  is  important  to  note  that  the  integration  in  Eqs.  (113)  and 
(1 14)  is  over  the  parameterized  rough  surfaee,  i.e.,  ds'  =  \n\dxdy' .  Consequently,  using 

Eq.  (110)  and  substituting  Eqs.  (1 15)-(1 19)  into  Eqs.  (113)  and  (114),  one  obtains  the 
following  expressions: 


and 


Q^'^'^'dk[dk[dxdy 

(122) 

Q'''"'dk[dk[dxdy  , 

(123) 

where  q  =  (f  -  k'  j  +  q^y  +  q^z ,  L'  =  n'  M' ,  and  N' 


Inlj'.  Without  further 


simplifieations,  no  analytieal  expression  exists  for  the  far-field  veetor  potentials,  L  (r ) 


and  N  (r) ,  given  in  Eqs.  (122)  and  (123).  This  is  beeause  the  integrands  in  Eqs.  (122) 

and  (123)  are  eomplieated  funetions  of  surfaee  height  and  surfaee  slopes;  namely,  h' , 

,  and  with  respeet  to  the  integrals  over  the  parameterized  rough  surfaee.  One 

9.  See  Appendix  B  for  more  details. 
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typically  simplifies  these  integrals  using  the  stationary-phase  (SP)  approximation  [22, 

35],  i.e., 

A(qr')  =  0  (124) 

OX  ay 

As  a  result,  the  relationships  found  in  Eq.  (120)  simplify,  such  that 

K'---  h,--—,  (125) 

^z 

and  in  turn,  L'  ~  L  and  N'  ~  N  in  Eqs.  (122)  and  (123).  Similar  to  the  PO 
approximation,  the  SP  approximation  physieally  dictates  that  refieetion  from  the  rough 
surfaee  is  locally  specular  and  excludes  all  loeal  diffraction  effects  [22,  35]. 

4,3  Scattered  field  cross-spectral  density  matrix 

The  analysis  presented  here  develops  closed-form  expressions  for  the  scattered 
field  CSDM  W*(rj,r2).  In  general,  W^(rj,r2)  depends  on  the  seattered  field  E^(r), 
sueh  that  in  the  far  field 

W'(r„r,)  =  (E'(r,)E‘t(r,)) 

_'(£;(r,)£r(r,))  (£;  (r,  )<  (r,  ))1 

=  l^pq(fi’*'2)  (p  =  6',^;  q  =  6',^) 


Using  Eq.  (112),  one  determines  the  matrix  elements  found  in  Eq.  (126)  as 

=  z  Z  [(i-*)(4'j)(A(ri)^U*'2)) 

j=x,y,z 

+  ^o(^'l)(02-j)(A(ri)^i(r2)) 

+  ^o(erl)(4'j)(^.(ri)^K*'2)) 

+  ^o(0i-i)(e2-j)(iV((rOiV;(r2))“ 


(127) 
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(128) 


^^;(ri,r2)  =  Qi2  X  Z  [-(^■*)(02-j){A(ri)^K*'2)) 

i=x,y,z  j=x,y,z 

+  ^o(i'j)(4'j)(A(ri)A^n*'2)) 

-%(erj)(02-j)(^.(ri)^’(*'2)) 

+  ^o(0rJ)(4'j)(^^i(ri)A^- (*■2)) 


fF;,(r„r,)  =  Q,,X  Z 

i=x,y,z  j=x,y,z 


-(01 ■])(!,  (r,)l*(r2)) 

-;7o(0ri)(e2-j)(A(ri)A^n*'2)) 

+  ^o(^'j)(4'j)(^.(ri)^K*'2)) 

+  ^o  (^■j)(02-j)(^.(ri)A'n*'2)) 


(129) 


and 


K('-r'-,)  =  ii„  Z  Z 

i=x,y,z  j=x,y,z 


(ei-i)(02-j)(A(ri)i;(r2)) 

-^o(6rj)(4-j)(A(ri)A^n*'2)) 

-^o(^'j)(02-j)(^'.(ri)^K*'2)) 

+  ^o(i'j)(4'j)(^'.(ri)A^n*'2)) 


(130) 


where  Q.^^=kl- — - .  (131) 

(4^) 


In  addition,  using  Eqs.  (1 13)-(125)'°,  one  determines  the  element-based  eorrelations 
found  in  Eq.  (130)  from  the  following  relationships: 


10.  In  using  Eqs.  (1 13)-(125),  one  must  assume  that  all  observation  points  are  in  the  far 
field  (ef  Eootnote  8,  p.  51),  the  physieal-opties  approximation  holds  (cf  Appendix  B), 
and  the  effeets  of  shadowing/masking  and  multiple  seattering  are  negligible  [38,  42]. 
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(L(r,)Lt(r,))  =  ^ 


jjunn  ,  (132) 

-L  -L  -L  -L  -oo  -<»  -oo  -<» 

Q-j'Kir,  ^jKir,  Qj‘lx\X\  Q~J‘ly2y2  IqJ‘Iz\K  g-y?z2*2  \ 


dk[^dk[2dk[^dk[2dx[dx2dy[dy2 


(L(r,)N’(r,)}  =  -— ^ 

{2n)  % 


iijinn  n;  _  (,33) 


-iKi^s  ^JKiKs  ^-Mx2^2  J^y\y\ 


dk[]  dk^  2  dkl  ^  dk^  ^  dx[dx\  dy\dy. 


(N(r,)L*(r,))=--^ 
{Itt)  j]. 


^r{X{K,.K,)T'{K,.K^))t\  ,  (134) 


-L  -L  -L  -L  - 


g-Xi'i  QjKyr,  ^NxiA  QJ'lyXl  I^qMxxK  'j 

jrJ  jiJ  jrJ  jrJ 


dk^^dk'^^dk^^dkl^dx^dx^dy^dy^ 


(N(r,)N’(r,))  = 


(24-)' 77.^ 


ijjinn  N,  (T'(i;„*;,)T"(<7;„^;j))-N;  ,  (bs) 


-JK\Ks  ^JK2Ks  ^-J<3x2^2 


g  J'-vl'5  qJ-v2-s  Qj-nxi-'i  g  J'ixl-'l  g  J^y2y2  j^j^zxK  g  J^z2^2  ^ 

dk[]  dk^  2  dkl  j  dk^  2  dx[dx'2  dy[dy. 


where  r;  =  .  Inherent  in  Eqs.  (132)-(135)  is  the  assumption  that  the  ineident  field 


plane-wave  speetrum  is  statistieally  independent  of  the  rough  surfaee.  This  assumption  is 
physieally  intuitive;  thus,  Eqs.  (132)-(135)  contain  two  separate  correlations. 
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The  first  correlation  is  with  respect  to  the  incident  field  plane-wave  spectrum. 

This  correlation  is  equivalent  to  a  dyadic  [cf.  Eq.  (Ill)];  namely, 

,  (136) 

where  O'  =  O'  ( A:]j ,  ,  k^^ ,  A:^2 )  •  The  second  correlation  is  with  respect  to  the 

parameterized  rough  surface.  This  correlation  is  a  joint  characteristic  function  of  the 

random  variables  h[  =  h  [x', y')  and  =  A  (xj , yQ ,  such  that 


(137) 


In  practice,  one  must  choose  a  form  for  this  joint  characteristic  function. 

A  very  common  choice  for  the  statistical  distribution  of  the  rough  surface  is  to 
assume  that  the  surface  heights  are  Gaussian  distributed  and  Gaussian  correlated.  In  so 
doing,  the  joint  probability  density  function  p' =  p[h[,h2)  of  the  random  variables  h[ 


and  h'^  takes  the  following  form  [36]; 


P  = 


exp 


2o-;(i-r'') 


(138) 


where  Y'  =  Y [x[  - x'^, y' - yQ  is  the  surface  autocorrelation  function,  such  that 


r'  =  exp 


1 

yT 

r 

K> 

1 

{yi-y'if 

1 

1 

1  _  . 

exp 

- 1 

<N 

1 

_ _ 1 

(139) 


History  shows  that  one  typically  chooses  Gaussian-Gaussian  (G-G)  models  for  analytical 
convenience  [93];  however,  other  models  exist  in  practice.  For  example,  the  stretched 
exponential-stretched  exponential  (SE-SE)  model  better  characterizes  surfaces  roughened 
by  random  industrial  processes  [138].  Basu  et  al.  highlighted  this  point  with  profilometer 
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measurements  of  sandblasted  metallie  surfaees  [139,  140].  Unfortunately  no  general 
analytieal  form  exists  for  the  SE  joint  oharaeteristie  funetion;  nevertheless,  the  analysis  of 
Basu  et  al.  also  showed  that  G-G  models  were  still  fairly  good  approximations  for 
sandblasted  metallie  surfaees  [139,  140].  Thus,  Fourier  transforming  the  joint  probability 
density  funetion  p'  in  Eq.  (138)  yields  the  desired  joint  oharaeteristie  funetion  for  the 
present  analysis  [36],  i.e.. 


z'=j  jp' 


:  exp 


(140) 


where  z'  =  z(Ki’K2’Ki’K2’  “  -^2  ’  “  U )  •  Note  that  throughout  the  literature, 

numerous  other  surfaoe  models  exist  in  addition  to  G-G  and  SE  models''. 

Using  the  relationships  found  in  Eqs.  (136)-(140),  the  integrands  in  Eqs.  (132)- 
(135)  still  oontain  oomplioated  funotions  with  respeot  to  the  souroe  and  surfaoe 
parameters.  To  simplify  the  analysis,  one  oan  separate  these  oomplioated  functions  into 
amplitude  and  phase  terms,  viz.. 


(L(r,)L^(u)) 


1 


L  L  L  L 


JJiJ 


dk^  1  2  1  2 


-L  -L  -L  -L  V-“ 


dxldx^dyldy^ 


(141) 


11.  For  example,  a  recent  publication  explored  the  use  of  non-Gaussian  surface 
autocorrelation  functions  [141]. 
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Here, 

(145) 

(146) 

f"=/(N, 

(147) 

and 

4>'  n;) 

(148) 

are  amplitude  dyadies  that  contains  all  of  the  amplitude  terms,  and 

g  =  [r^  +  (y .  v) _y; ] - [r,  +  (y  ■  v) 

+ -  y'lKi )  ( i49) 

Kq  Kq 

+  (x-fX-x-f2x;)  +  (y-fiy;-y-f2y;) 


is  a  common  phase  function  that  contains  all  of  the  phase  terms.  Without  further 
simplifications,  no  closed-form  expressions  exist  for  the  integral  relationships  given  in 
Eqs.  (141H144). 
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To  simplify  the  integrals  found  in  parenthesis  in  Eqs.  (141)-(144),  the  analysis 
uses  an  asymptotie  mathematical  technique  known  as  the  method  of  stationary  phase 
(MoSP)  [26,  27,  37].  In  so  doing,  one  assumes  that  the  amplitude  terms  are  slowly 
varying  in  the  interval  (-00,00) .  One  must  also  assume  that  the  phase  terms  are  rapidly 

oscillating  in  the  interval  (-00,00)  except  near  special  points  where  the  rate  of  change  is 

zero  or  “stationary.”  These  special  points  are  called  critical  points  of  the  first  kind  [26]. 
Away  from  these  points,  the  phase  terms  are  rapidly  oscillating  and  the  positive  and 
negative  contributions  of  the  integrand  in  Eq.  (144)  effectively  cancel  out. 

Using  the  MoSP  to  simplify  Eqs.  (141)-(144)  has  two  implications  with  regards  to 
the  macro-scale  scattering  geometry  given  in  Eigure  9a.  The  first  implication  is  with 
respect  to  the  v  component  of  the  incident  propagation  vector  k' ;  namely, 

(150) 

In  particular,  the  analysis  assumes  that  »  k[  and  ;  as  a  result. 


ing 


(151) 


This  physically  implies  that  the  incident  electromagnetic  fields  are  highly  directional 
being  predominately  directed  along  the  v  direction  in  Eigure  9a.  The  second  implication 
is  that  the  distance  from  the  source-plane  origin  to  the  surface-plane  origin  must  be  much, 
much  greater  than  half  the  surface  length,  i.e.,  »  E  ,  which  is  typically  the  case  for 


most  directed-energy  and  remote-sensing  engagement  scenarios.  To  provide  some  idea 
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of  how  much  greater,  letting  r^=5L  ,  r^=  25 L  ,  and  =  lOOZ  results  in  pereentage 

errors  of  18%  ,  4%  ,  and  1%  ,  respeetively.  As  such,  using  the  MoSP  to  simplify  Eqs. 
(141)-(144)  results  in  the  following  relationships: 


(l  (e  )L^'  (rj ))  ~  — — — C , 

(152) 

(L(r,)N'(r2))-  °  ^ 

[Itt)  //ofi 

(153) 

(154) 

and 

(155) 

where  C ,  K,  A4 ,  and  A^  are  dyadics  that  eontain  all  of  the  amplitude  and  phase  terms 

12 

evaluated  at  the  eritieal  points  of  the  first  kind  which  one  determines  as 


^o(y-«) 


^uL2 


■yi2 


(156) 


The  analysis  explieitly  defines  C  ,  K ,  A4 ,  and  M  in  Appendix  D  for  different  material 
substrates,  i.e.,  dieleetries,  eonduetors,  and  a  perfeet  eleetrieal  eonduetor  (PEC). 

Provided  Eqs.  (152)-(156)  and  Appendix  D,  one  is  still  left  with  integrals  with 
respeet  to  the  parameterized  rough  surfaee.  These  integrals  take  the  following  element- 
based  form: 


12.  See  Appendix  C  for  more  details. 
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L  L  L  L 


mn  J  J  J  J 

-L  -L  -L  -L 


K  '  K  '  ^o(y-u)  /  ^o(y-u)  / 

— Xj,— Xj, - y„ - J2 


v'; 


K  '  K  >  ^o(y-«)  /  ^o(y-»)  ,  ,  ,  ,  , 

Xp  — Xj, - Ji, - yy,x^-x^,y^-y^ 

r,  r 


v':. 


exp [7I0  (x  ■  yx[  -  X  ■  r^x'^ )]  exp [jk^  (y  ■  r,y[  - y  ■  1^272 )]  . 
exp 


(157) 


1 

(N 

1 _ 

.^o(y'U)%  /2  /2^ 

exp 

J  2^^  ) 

exp  [-7I0  (y  ■  v)  (X  -  y; )]  dx[dx^dy[dy[ 

where  (m  =  x,m;  n  =  x,m)  .  In  Eq.  (157),  is  equivalent  to  the  Fourier  transform  of 
the  incident  field  CSDM  elements  found  in  Eq.  (108),  i.e., 

<.=1111  WL  (Ppp2)e^'''''‘  dx,dx,du,du. 


_  ^'AnABmn 

(aL-bL) 


exp 


exp 


-d  {k\\  —d  {k\\  +2b  k\k\ 

mn  \  xl  J  mn  \  xl  J  mn  xl  x2 

-d  {k\\  —d  {k\\  +2b  k\k\ 

mn  y  u\  J  mn  y  u2  J  mn  u\  u2 


(158) 


where  (m  =  x,M;  n  =  x,M)  and 


^  u  u  1 


^mn  A  .  2  ‘  "^mn  "^mn 


4w( 


2i 


2  «mn=- 


=■ 


(159) 


4(a"  -b^  )  ™  4(a"  -b^  )' 

One  can  reduce  the  integrals  found  in  Eq.  (157)  into  closed-form  expressions.  For 
this  purpose,  the  analysis  first  performs  the  following  variable  transformations: 

Xj  =  x[-X2  x^=xi'+x2  y,=y'i-y2  ya=y'  +  /2,  (160) 

so  that  Eq.  (157)  simplifies  into  the  following  expression: 
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exp 


2L  2i-|>'„|  2L  2i-|x„ 


^K,+4) 


-2L  \ya\-2L  -2L  \x„\-2L 

kl 


J  J  j  J  exp 

Xa\~2L 

2^,  (<1™+*™) 


^o^^/,^zi^z2exp 


7  2 

f  2  M 

'  _ZfL 

exp 

V  J 

V  )_ 

exp 


exp 


exp 

exp 

exp 


^0  (y-uf 

2r" 


2r' 


exp 


-id  —b  ]x^ 

Y  mn  mn  J  a 


^0  (y-uf 

2r} 


id  —b  \y^ 

\  mn  mn  j  y  a 


7^(4 +  ^x2)^^ 


exp 


j-^iAx-^x2)^a 


j^(^yx  +  '^y2)yd  exp  7 y  ( ) 7. 


•  K 

1 

<N 

1 _ _ 

exp 

L  2r,  J 

exp 

7  _  yayd 

2r, 

dxjdxjy^dy^ 


(161) 


where  (m  =  x,M;  n  =  x,u)  and 

z^xi,2  =x-ri,2  z^2,i,2  =y-ri,2-y-v  4^2=z-ri,2-Z'V.  (162) 

From  here,  one  must  handle  the  exponential  term  eontaining  the  surfaee  autoeorrelation 
funetion  (i.e.,  the  first  exponential  term  inside  the  integrals  above).  History  shows  that 
there  are  two  separate  ways  to  go  about  this.  The  first  is  to  expand  the  said  exponential 
term  in  a  Taylor  series  and  proeeed  with  the  evaluation  of  the  integrals  [36]. 
Mathematieally,  this  approaeh  is  applieable  to  all  surfaees;  however,  beeause  the  series  is 
slowly  convergent,  the  analysis  limits  this  approach  to  smooth-to-moderately  rough 
surfaces — the  next  sub-section  develops  a  closed-form  expression  for  this  case.  The 
other  approach  involves  expanding  the  surface  autocorrelation  function  [cf.  Eq.  (139)]  in 
a  Taylor  series  and  retaining  only  the  first  and  second  order  terms  [36].  This  treatment  is 
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applicable  to  very  rough  surfaees — sub-section  4.3.2  develops  a  closed-form  expression 
for  this  case. 


4.3.1  Smooth-to-moderately  rough  surfaces 

When  considering  smooth-to-moderately  rough  surfaces,  one  must  expand  the 
exponential  term  containing  the  surface  autocorrelation  function  found  in  Eq.  (161)  in  a 
Taylor  series.  Specifically, 


exp 

kWhAA2  exp 

r 

j)2 

exp 

1 

z' 

^  -a  u 
^1  ^ 

K  ^h) 

E- 

z=o 


/! 


exp 

(  1  ^ 
^  2 

-fT^d 

exp 

^  ^ 
‘  ^  2 

1  ^2  y^ 

K  H  J 

K  ) 

(163) 


Substitution  of  Eq.  (163)  into  Eq.  (161)  allows  one  to  then  separate  the  integrals  over  the 
parameterized  rough  surface,  i.e.. 


mn  A  (  1  j  2  \  ^ 

4  a  -b  I 

\  mn  mn  J 


Z- 

/=0 


/! 


21 

J 


exp 


-2L 

2i-|j 


-id  —b  ) 

\  mn  mn  j 


I  exp 


,\-2L  V 

f 

exp 


2r; 


^  ^ 
<  2 

•h  J 


exp 


j^iAx-A2Va 


exp 


2r' 


[a  +b 

mn  mn  J 


(^.1  +  ^x2)] 

V  J 


dx^dx^ 


2L 

f 

^0  y-«  i~  1  \  2 

- 

J  ^^P 

-2L 

.  2  («mn  Kn)yu 

2f  ^  ^  J 

exp 

j^[^yl-'^y2)ya 


2Z,-|v^| 

f 

|>’»|-2i  V 
( 


^  I  ^ 
^  2 

J 


exp 


^0  (y-»)' 
'IrJ 


(«mn+^mn) 


2 

yd 


exp 


-J 


'^0 

2r 


(y-^fya-^s(^yl+^y2) 


yd 


dyddy. 


(164) 
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This  is  a  very  important  step  in  the  analysis.  It  allows  for  the  development  of  a  elosed 
form  expression  for  smooth-to-moderately  rough  surfaee  eonditions  without  having  to 
convert  to  polar  coordinates. 

It  is  a  relatively  straight  forward  process  to  evaluate  the  integrals  over  and 

in  Eq.  (164);  however,  complex  error  functions  result  due  to  the  parameterized  rough 
surface  [10,  11].  These  complex  error  functions  model  diffraction  caused  by  the  incident 
radiation  over-illuminating  the  rough  surface  and  are  negligible  under  certain  conditions. 
As  such,  the  present  goal  in  the  analysis  is  to  determine  conditions  in  which  one  can 
extend  the  and  limits  of  integration  to  (-00,00)  in  Eq.  (164).  The  necessary 


conditions  for  these  approximations  occurs  when 


exp  --^x]  exp 


/  2 

exp  -^yd  exp 
K  J 


^0 


id  +b  \yl  >  S 

Y  mn  mn  J  y  a  y 


(165) 


where  and  are  user-defined  parameters  and  denote  the  points  at  which  the 
exponential  functions  with  respect  to  x^  and  y^  no  longer  maintain  significant  value  in 
Eq.  (164).  Eurthermore,  if  /  =  0  and  x^  =  y^  =  2L  ,  the  arguments  of  the  exponentials  are 
at  a  minimum  and  the  x^  and  y^  limits  of  integration  are  at  a  maximum.  One  then 


derives  the  following  conditions  from  Eq.  (165): 


(166) 
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These  eonditions  physieally  mean  that  the  projeeted  fully  eoherent  ineident  beam  size 
must  “fit”  on  the  rough  surfaee.  In  satisfying  these  eonditions  the  and  limits  of 


integration  extend  to  (-00,00)  in  Eq.  (164).  Note  that  if  =  S ,  the  seeond 


eondition  in  Eq.  (166)  beeomes  more  stringent  sinee  y  •  u  =  -eos(^')  <1  [ef  Eigure 


9a].  One  then  determines  how  well  the  ineident  beam  “fits”  by  S — the  smaller  the  S , 
the  more  aeeurate  the  approximation.  Also  note  that  if  the  projeeted  fully  eoherent 
ineident  beam  size  does  not  “fif’  on  the  rough  surfaee,  eomplex  error  funetions  result  [10, 
11],  and  one  has  to  evaluate  the  follow-on  integral  expression  numerieally. 

Assuming  that  the  eonditions  in  Eq.  (166)  hold,  subsequent  evaluation  of  the 

1 3 

integrals  over  and  simplifies  Eq.  (164),  so  that 


13.  One  must  eomplete  the  square  in  the  exponential  terms  and  use  of  the  following 
integral  relationship  [125,  p.  266];  |  Qxp(^-at^^Qxp[-jbt)dt  =  exp[^-6Y(4a)J  , 

where  a  >  0  . 
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exp 
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where  (m  =  x,w;  n  =  x,u)  and 


2?™=^i(’U«.„+6™)(y  flf+2fr; 


(168) 


Similar  to  the  analysis  presented  above,  one  ean  extend  the  x^  and  limits  of 
integration  to  (-00,00)  in  Eq.  (164)  and  subsequently  Eq.  (167).  Here,  the  neeessary 


conditions  for  these  approximations  occurs  when 
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(169) 


where 


^mn 


(170) 


is  the  element-based  source  ratio,  and  and  S  are  again  user-defined  parameters. 


respectively.  They  denote  the  points  at  which  the  exponential  functions  with  respect  to 
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and  no  longer  maintain  significant  value  in  Eq.  (167).  The  conditions  given  in  Eq. 
(169)  physically  mean  that  the  projected  partially  coherent  incident  beam  size  must  “fit” 
on  the  rough  surface.  In  satisfying  these  conditions  the  and  limits  of  integration 

extend  to  (-00,00)  in  Eq.  (167).  Additionally,  if  5^  =  5^  =  5 ,  the  second  condition  in  Eq. 

(169)  becomes  the  most  stringent  within  the  smooth-to-moderately  rough  surface 
analysis. 

Assuming  that  the  conditions  in  Eq.  (169)  hold,  one  can  then  evaluate  the 
remaining  integrals  in  Eq.  (167)'"^.  In  so  doing,  the  following  closed- form  expression 
results: 
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where  (m  =  x,M;  n  =  x,M)  and 


14.  See  Eootnote  13,  p.66. 
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(172) 


AL=rs^\  (y-uf  -Kn)^L 

At  first  glance,  the  summation  somewhat  obseures  the  physieal  interpretation  of  Eq. 
(171).  However,  further  examination  of  this  elosed-form  expression  shows  that  the 
exponential  terms  on  the  third  and  fourth  lines  generally  drive  the  angular  extent  of  the 
speetral  density  (SD),  whereas  the  exponential  terms  on  the  fourth  and  fifth  lines 
generally  drive  the  angular  extent  of  the  speetral  degree  of  eoherenee  (SDoC).  Before 
exploring  these  points  further  in  the  next  ehapter,  the  analysis  considers  very  rough 
surfaees'^. 


4.3.2  Very  rough  surfaces 

When  eonsidering  very  rough  surfaees  in  the  analysis,  one  must  expand  the 
surfaee  autoeorrelation  funetion  found  inside  the  first  exponential  term  in  Eq.  (161). 
Here,  one  retains  only  the  first  and  seeond  order  terms  in  a  Taylor  series.  To  make  this 
eoneept  manifest,  the  analysis  first  writes  the  joint  eharaeteristie  funetion  found  in  Eq. 
(140)  in  an  alternative  form,  where 
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In  the  ease  of  very  rough  surfaee  eonditions,  i.e.. 


A- 
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(173) 


15  The  following  eriterion;  >  O.S/Iq  ,  helps  in  diseeming  the  transition  point  from  the 

smooth-to-moderately  rough  surfaee  regime  to  the  very  rough  surfaee  regime  and  is  an 
empirieally  determined  relationship  within  the  analysis. 
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kWlAAi » 1 . 


(174) 


the  alternative  form  found  in  Eq.  (173)  maintains  signifieant  value  when 
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Sinee  all  of  the  observation  points  of  interest  are  in  the  far  field,  if  one  then  eonsiders  that 


z7., 


(176) 


2^5»2  24 

the  relationship  found  in  Eq.  (175)  is  only  possible  for  small  and  . 

With  Eqs.  (173)-(176)  in  mind,  it  makes  sinee  to  expand  the  exponential  funetions 
found  in  Eq.  (175)  and  retain  only  the  first  and  seeond  order  terms,  so  that 
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(177) 


Substituting  Eq.  (177)  into  Eq.  (161)  allows  one  to  again  separate  the  integrals  over  the 
parameterized  rough  surfaee.  As  mentioned  before,  this  allows  for  the  development  of  a 
elosed-form  expression  for  very  rough  surfaee  eonditions  without  having  to  eonvert  to 
polar  eoordinates.  Carrying  out  the  subsequent  integrations'^,  the  following  elosed-form 
expression  results: 


16.  See  Eootnote  13,  p.  66. 
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where  [m  =  x,u;  n  =  x,u)  and 


=iX  (S^  +L}  +  2tkA;A,i^,2 

^  I’..  I’.: 

-C  =A<i+4(a^,-6^)2r^ 


4;.='f(;(y")  +4(5™-*™)^’.: 


(179) 


In  using  the  elosed-form  expression  given  in  Eq.  (178),  the  analysis  must  satisfy  the 
eonditions  found  in  Eq.  (169). 

The  elosed-form  expression  obtained  in  Eq.  (178)  is  remarkably  physieal.  Eor 
instanee,  the  exponential  terms  on  the  seeond  and  third  lines  of  Eq.  (178)  are 
predominately  responsible  for  the  angular  extent  of  the  seattered  SD.  These  exponential 
terms  are  funetions  of  the  sum  of  the  squares  of  the  observation  projeetions,  i.e.,  2  and 

2 .  On  the  other  hand,  the  exponential  terms  on  the  fourth  and  fifth  lines  of  Eq.  (178) 
determine  the  angular  extent  of  the  seattered  SDoC.  Note  that  these  terms  are  funetions 
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of  the  difference  of  the  observation  projections,  i.e.,  “^v2|  •  Thus,  one 

can  state  that  the  closed- form  expression  obtained  in  Eq.  (178)  allows  the  scattered  field 
CSDM  to  maintain  its  GSM  form  with  respect  to  2  and  z?  j  2 .  The  analysis  examines 

these  points  more  closely  in  the  next  chapter. 
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5  Exploration  of  the  3D  vector  solution 


The  purpose  of  this  ehapter  is  to  explore  the  3D  veetor  solution  obtained  above  in 
Chapter  4.  As  stated  in  Chapter  1,  by  formulating  the  3D  veetor  solution  in  a  manner 
eonsistent  with  Wolfs  unified  theory  of  coherenee  and  polarization  [12,  13],  all  physieal 
implieations  inherent  in  Wolfs  work  apply  here.  Aeeordingly,  one  ean  readily  formulate 
the  seattered  speetral  degree  of  eoherenee  (SDoC)  (ri,r2 )  >  the  normalized  seattered 


speetral  density  (SD)  5*^  (r) ,  and  the  seattered  degree  of  polarization  (DoP)  C*'  (r)  from 
the  elosed-form  expressions  developed  above  for  the  seattered  field  eross  speetral  density 
matrix  (CSDM)  (fi  A2)  •  The  analysis  uses  the  following  relationships  [13]; 


Tr{w^(fi,r,)} 


Tr{w^(r„rO}JTr{w^(r„r,)} 


(180) 


Tr{w"(r,r)} 
max  Tr|w*(r,r)j 


(181) 


and 


F'(r)  = 


4Det{w"(r,r)} 
(Tr{w^(r,r)})'  ’ 


(182) 


where  Tr{---|  denotes  the  traee  operation,  Det{---|  denotes  the  determinant  operation. 


and  r  =  2  corresponds  with  a  single  observation  point.  These  relationships  eontain 

measurable  quantities  in  praetiee  and  serve  as  metries  in  whieh  to  eompare  the  3D  veetor 
solution  to  previously  validated  solutions  and  empirieal  measurements. 

Mueh  of  the  analysis  presented  in  this  chapter  uses  a  5.08  cm  x  5.08  cm 
Labsphere  Infragold  coupon  [142].  It  also  uses  a  nominal  far-field  setup,  where 
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/Iq  =  1 0.6  |am ,  =  Tj  j  =  1 85  cm ,  and  w^=\.9  mm  .  As  such,  the  Labsphere  Infragold 

coupon  maintains  the  following  complex  index  of  refraction;  n  =  13.45-763.62  [143]. 
Note  that  a  KLA  Tencor  Alpha-Step  IQ  Surface  Profiler  [144]  determined  the  surface 
statistics  of  the  Labsphere  Infragold  coupon  as  <7^  =  1 1 .09  pm ,  f  ^  =  1 16.9  pm ,  and 
(7^.  =  0.2441  rad  [cf.  Eq.  (191)]  using  four  1  cm  scans  (step  size  0.2  pm).  These  surface 

statistics  relate  to  very  rough  surfaee  conditions  [cf.  Eq.  (174)]. 

5.1  Comparison  with  the  2D  scalar-equivalent  solution 

In  order  to  eompare  the  3D  veetor  solution  to  the  previously  validated  2D  sealar- 
equivalent  solution  [10,  11],  the  analysis  assumes  horizontally  polarized  (s-pol) 
illumination  and  an  in-plane  scattering  geometry,  i.e.,  =  0  and 

(j)'  =  =  90  .  This  provides  the  setup  needed  to  make  a  fair  comparison  between  the 

two  solutions. 

5.1.1  Angular  spectral  degree  of  coherence  radius 

An  individual  can  formulate  a  closed-form  expression  that  describes  the  angular 
extent  over  whieh  the  scatter  field  is  correlated,  i.e.,  the  angular  SDoC  radius.  In  general, 
the  angular  SDoC  radius  provides  a  gauge  for  the  average  speckle  size  observed  in  the  far 
field  and  is  a  quantity  of  importance  when  dealing  with  direeted-energy  and  remote¬ 
sensing  applications.  Note  that  the  analysis  presented  here  is  highly  analogous  to  that 
performed  for  the  2D  scalar-equivalent  solution  [10,  11]. 

Because  of  the  summation  in  Eq.  (171),  it  is  not  possible  to  derive  a  closed-form 
expression  for  the  angular  SDoC  radius  for  smooth-to-moderately  rough  surfaces.  Thus, 
the  present  analysis  is  limited  to  very  rough  surfaces.  Assuming  that  Eq.  (169)  holds,  so 
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that  the  incident  illumination  “fits”  on  the  rough  surface,  the  exponential  terms  on  the 
fourth  and  fifth  lines  of  Eq.  (178),  in  general,  determine  the  angular  extent  of  the 
scattered  SDoC.  Provided  s-pol  illumination  and  an  in-plane  scattering  geometry,  only 
exists  within  the  analysis,  and  the  difference  of  the  observation  projections  simplify, 
such  that 


Ax-Ai  = 


sm 


[O^ )  cos  )  -  sin  (^2^ )  cos  {(pl )  =  0 


(183) 


and 


sm 


(^j')sin(^*)-sin(^2  )sin(^2 )  =  sin^^,*  )-sin(^2 )  •  (184) 


Consequently,  the  following  “correlation”  exponential  y  results  from  Eq.  (178)  for  : 


Y=  exp 


exp 


2„2 


(y-fi)  4; 


S  XX 


*, )  [sin  ( s; )  -  sin  («; ) 


2(y'ii)  4; 


sm 


(<?;)-sin(e‘)T 


(185) 


Upon  setting  y  equal  to  1/e,  the  following  expression  results: 


sin(^1-sin(^j|^^ 
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Because  the  magnitude  of  the  argument  of  y  is  large  (specifically  the  klr^ 
term),  sin [O^ )  ~  sin(^2 )  Eq.  (185)  to  have  a  significant  value.  This  implies  that 

01  ~  61  and  that  y  is  approximately  a  function  of  A^'  =  61  -6[  .  Using  this  insight,  Eq. 
(186)  simplifies  because 
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(187) 


sin  ( 0^ )  =  sin  ( ^2  +  ^  ) 

=  sin  (^2' )  cos  ( A^’ )  +  cos  (^2 )  > 

~  sin  (^2 )  +  cos  {Ol )  A^' 

cos  ( )  =  cos  ( ^2  +  ) 

=  cos(^2)cos(A^')-sin(^2  )sin(A0’) ,  (188) 
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After  some  simple  algebra,  the  expression  for  the  angular  SDoC  radius  beeomes 

1 


A^-' 


ll/e 


8w,^  ^  2 

l  +  (2/«„f 

\  !  XX/  2  2/1 

^oCr,.(l  +  57) 

r 

1  + 

i 

^2^ 

2^ 

_rT,^,(l  +  57)_ 

(190) 


where  =^xc/^^  Is  ^  souree  ratio  [of  Eq.  (170)],  (D  =  oos(^2)/cos(^')  is  aprojeotion 
ratio,  =  wjr^  is  the  souree  half  angle  (viewed  from  the  rough  surfaoe),  and 


is  the  surfaoe  slope  standard  deviation  [36]'^. 


(191) 


17.  Based  on  the  assumptions  used  within  the  analysis  [of  Footnote  10,  p.56],  valid 
surfaoe  slope  standard  deviations  must  satisfy  the  following  oondition;  <  0.25  rad 
[35,  66]. 
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For  all  intents  and  purposes,  one  can  neglect  the  term  involving  in  Eq. 


(190).  This  ratio  results  in  values  on  the  order  of  10  for  most  directed-energy  and 
remote-sensing  engagement  scenarios.  With  this  said,  one  can  also  claim  that  the  source 
term  contained  in  the  radical  above  is  much  greater  than  the  surface  term.  Thus, 
factoring  out  the  source  term  and  using  the  binomial  approximation  yields 
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I  ll/e  ^ 
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(192) 


For  most  cases  of  interest,  one  can  neglect  the  second  term  contained  within  the 
parenthesis  in  Eq.  (192).  It  only  provides  a  “small”  correction  to  the  angular  SDoC 
radius  due  to  the  surface  parameters.  As  a  result,  the  angular  SDoC  radius  becomes  a 
function  of  only  the  source  parameters.  This  is  highly  analogous  to  the  result  obtained  by 
the  2D  scalar-equivalent  solution  [10,  1 1].  It  is  also  consistent  with  the  classic,  narrow- 
band,  fully  coherent  illumination  result  derived  by  Goodman  [90]. 

S.I.2  Angular  spectral  density  radius 

An  individual  can  also  formulate  a  closed-form  expression  for  the  angular  SD 
radius.  In  general,  the  angular  SD  radius  provides  a  gauge  for  the  size  of  the  average 
power  distribution  observed  in  the  far  field.  This  is  a  quantity  of  importance  when 
dealing  with  directed-energy  and  remote-sensing  applications.  The  analysis,  yet  again,  is 
highly  analogous  to  that  performed  for  the  2D  scalar-equivalent  solution  [10,  11]. 
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As  was  the  case  for  the  angular  SDoC  radius,  the  summation  in  Eq.  (171)  does  not 
allow  for  a  closed-form  expression  for  the  angular  SD  radius  for  smooth-to-moderately 
rough  surfaces.  Furthermore,  one  must  limit  the  analysis  to  near  normal  incidence,  so 

that  |y  ■  u|  =  cos  )  ~  ^  •  Assuming  that  Eq.  (169)  holds,  so  that  the  incident  illumination 

“fits”  on  the  rough  surface,  the  exponential  terms  on  the  second  and  third  lines  of  Eq. 
(178)  predominantly  determine  the  angular  extent  of  the  scattered  SD  for  very  rough 
surfaces.  Provided  s-pol  illumination,  only  exists  within  the  analysis.  In  addition, 

for  an  in-plane  scattering  geometry  and  a  single  observation  point,  i.e.,  r  =  Tj  j ,  the 
observation  projections  simplify,  such  that 


A  =  =  sin  )  cos  [f )  =  0 

(193) 

and 

dy  =  dy^  2  =  sin  ) sin  )  =  sin  ) . 

(194) 

Consequently,  the  following  “power-distribution”  exponential  (3  results  from  Eq.  (178) 
for  : 


/5  =  exp 


mn  Tnn 


(195) 


Upon  setting  f5  equal  to  1/e ,  the  following  expression  results: 


l-l-cos(^^Jl  -I--  ^ 


Iklwl 


(196) 


Since  sin^  ^  (^i/e )  ’  (196)  manipulates  into  a  quadratic  equation  in  terms 


of  cos(^/g) ,  where 
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Solving  this  quadratic  equation  (only  the  positive  root  makes  physieal  sense),  the  angular 
SD  radius  beeomes 


For  most  oases  of  interest,  the  souroe  and  surfaoe  terms  oontained  within  the  radioal  in 
Eq.  (198)  are  negligible.  They  only  provide  a  “small”  oorreotion  to  the  angular  SD  radius 
due  to  the  souroe  parameters.  Thus,  the  angular  SD  radius  beeomes  a  funotion  of  only  the 
surfaoe  parameters.  This  is  highly  analogous  to  the  result  obtained  by  the  2D  soalar- 
equivalent  solution  [10,  11]. 

5.1.3  Fully  coherent  illumination  validation 

In  order  to  validate  the  angular  radii  developed  above,  the  present  analysis  uses 
the  Labsphere  Infragold  ooupon  and  the  nominal  far-field  setup  (desoribed  above)  with 
fully  ooherent  illumination  at  normal  inoidenoe,  so  that  =  lOOOw^  and  O'  =0° . 


Provided  this  setup.  Figure  10  shows  a  comparison  between  the  3D  veotor  solution,  the 
2D  soalar-equivalent  solution  [10,  11],  and  a  full-wave  2D  method  of  moments  (MoM) 
solution  [139,  140].  The  2D  MoM  solution  obtained  the  seattered  field  from  400 
independent  rough  surfaoe  realizations  simulated  using  the  method  desoribed  by  Yura 
and  Hanson  [145]  with  a  Gaussian-Gaussian  (G-G)  probability  distribution  funotion 
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(PDF)  [cf.  Eq.  (138)]^^.  Note  that  the  results  mateh  up  well  and  that  the  elosed-form 
expressions  for  the  angular  radii  behave  as  predicted.  For  instance,  the  dashed  vertical 
line  in  Figure  10a,  which  represents  the  angular  SDoC  radius  [cf  Eq.  (192)],  identifies 
the  correct  1/e  location,  whereas  the  dashed  vertical  lines  in  Figure  10b,  which  represent 


the  angular  SD  radius  [cf  Eq.  (198)],  come  close  to  the  1/e  locations.  The  analysis 
performed  in  the  next  sub-section  further  explains  this  small  disagreement  with  respect  to 
the  angular  SD  radius;  nonetheless,  the  results  in  Figure  10  help  to  validate  the  3D  vector 
solution  assuming  fully  coherent  illumination  and  very  rough  surface  conditions. 


Figure  10.  Comparison  between  a  full-wave  2D  method  of  moments  (MoM) 
solution,  the  2D  scalar-equivalent  solution,  and  the  3D  vector  solution  for  fully 
coherent  illumination  at  normal  incidence  of  a  very  rough  conducting  surface,  (a) 
shows  the  magnitude  of  the  scattered  spectral  degree  of  coherence  as  a  function  of 
the  difference  between  two  polar  angles,  whereas  (b)  shows  the  normalized  scattered 
spectral  density  as  a  function  of  a  single  polar  angle. 


18.  A  full-wave  3D  MoM  solution  is  unrealizable  at  optical  wavelengths  with  the  current 
fully  coherent  setup — the  computational  sampling  and  memory  requirements  are  far  too 
great  for  modem  desktop  computers. 
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In  order  to  validate  the  3D  vector  solution  for  smooth-to-moderately  rough 
surface  conditions,  the  analysis  assumes  the  same  nominal  far-field  setup  used  above  but 
varies  the  surface  height  standard  deviations  of  the  Labsphere  Infragold  coupon,  so  that 
<7^  =  0.0\,  O.IAq,  and  0.3Aq  .  Figure  1 1  shows  the  results  for  this  setup.  Note  that  for 
both  the  magnitude  of  the  scattered  SDoC  (Figure  11a)  and  the  normalized  scattered  SD 
(Figure  1  lb),  the  results  deviate  for  small  <7;,  between  the  3D  and  2D  solutions.  This  is 

most  likely  due  to  fact  that  the  2D  scalar-equivalent  solution  is  confined  to  a  single  plane, 
whereas  the  3D  vector  solution  is  not.  The  results  of  Hyde  et  al.  showed  excellent 
agreement  between  the  2D  scalar-equivalent  solution  and  a  full-wave  2D  MoM  solution 
using  similar  setup  parameters  [11].  With  this  said,  the  2D  solutions  do  not  capture  all 
the  physics  related  to  the  3D  vector  problem — the  next  sub-section  examines  this  point 
further  at  normal  and  non-normal  incidence  using  partially  coherent  illumination. 


Figure  11.  Comparison  between  the  2D  scalar-equivalent  and  3D  vector  solutions 
for  fully  coherent  illumination  at  normal  incidence  of  smooth-to-moderately  rough 
conducting  surfaces,  (a)  shows  the  magnitude  of  the  scattered  spectral  degree  of 
coherence  as  a  function  of  the  difference  between  two  polar  angles,  (b)  shows  the 
normalized  scattered  spectral  density  as  a  function  of  a  single  polar  angle. 
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5.1.4  Partially  coherent  illumination  validation 

The  present  analysis  varies  the  coherence  of  the  incident  illumination,  so  that 
a^=2,  1,  0.5,  and  0.25  [cf  Eq.  (170)],  where  w^=\.9  mm  and  =  oc^w^ .  These 

values  relate  to  a  coherent  source  and  to  a  relatively  incoherent  source,  respectively.  In 
addition,  the  present  analysis  varies  the  surface  roughness  of  the  Labsphere  Infragold 
coupon,  so  that  =  0.01  rad ,  0.05  rad ,  0.1  rad ,  and  0.25  rad  [cf  Eq.  (191)],  where 

f  ^  =  1 16.9  pm  and  <7^  =  .  These  values  relate  with  smooth- to-very  rough 

surface  conditions,  respectively.  Eigure  12  and  Eigure  13  below  show  results  for  the 
magnitude  of  the  scattered  SDoC  and  the  normalized  scattered  SD  for  partially  coherent 
illumination  at  normal  incidence,  i.e.,  ff  =  Q° .  Note  the  excellent  agreement  between  the 
2D  and  3D  solutions  (similar  to  that  achieved  for  the  fully  coherent  illumination 
analysis). 

Eigure  14  shows  the  results  for  the  normalized  scattered  SD  for  partially  coherent 
illumination  at  non-normal  incidence,  i.e.,  6'  =  40° .  Here,  the  results  deviate  between  the 
2D  and  3D  solutions  as  (T;,.  increases.  This  deviation  is  deterministic  in  practice.  Eor  s- 

pol  illumination  and  an  in-plane  scattering  geometry,  the  2D  and  3D  solutions  have  the 
same  functional  dependence  in  the  exponential  term  which  predominately  drives  the 
angular  extent  of  the  scattered  SD.  The  analysis  explored  this  functional  dependence 
above  in  the  derivation  of  the  angular  SD  radius.  With  that  said,  there  is  a  small 
difference  contained  in  the  amplitude  terms  when  comparing  the  2D  and  3D  solutions. 
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Figure  12.  Comparison  between  the  2D  solution  (circles)  and  3D  solution  (lines)  for 
partially  coherent  illumination  at  normal  incidence  of  smooth-to-very  rough 
conducting  surfaces,  (a)-(d)  shows  the  magnitude  of  the  scattered  spectral  degree  of 
coherence  as  a  function  of  the  difference  between  two  polar  angles  for  varying 
source  parameter  ratios  and  surface  slope  standard  deviations. 
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Figure  13.  Comparison  between  the  2D  solution  (circles)  and  3D  solution  (lines)  for 
partially  coherent  illumination  at  normal  incidence  of  smooth-to-very  rough 
conducting  surfaces,  (a)-(d)  shows  the  normalized  scattered  spectral  density  as  a 
function  of  a  single  polar  angle  for  varying  source  parameter  ratios  and  surface 
slope  standard  deviations. 
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Figure  14.  Comparison  between  the  2D  solution  (circles)  and  3D  solution  (lines)  for 
partially  coherent  illumination  at  non-normal  incidence  of  smooth-to-very  rough 
conducting  surfaces,  (a)-(d)  shows  the  normalized  scattered  spectral  density  as  a 
function  of  a  single  polar  angle  for  varying  source  parameter  ratios  and  surface 
slope  standard  deviations. 


85 


The  3D  vector  solution  contains  a  ij amplitude  factor;  whereas,  the  2D 

scalar-equivalent  solution  contains  only  a  ij amplitude  factor  [cf  Eqs.  (172)  and 

(179)].  These  amplitude  factors  (in  addition  to  polarization  terms)  appear  in  front  of  the 
exponential  term  which  predominantly  drives  the  angular  extent  of  the  scattered  SD. 
Functionally,  these  amplitude  factors  tend  to  “push”  the  scattered  SD  to  the  right, 
whereas  the  polarization  terms  tend  to  “pull”  the  scattered  SD  to  the  left  (for  positive 
incident  angles — the  opposite  is  true  for  negative  incident  angles).  The  additional 
amplitude  factor  contained  in  the  3D  solution  pushes  the  scattered  SD,  so  that  the  peak  of 
the  far-field  power  distribution  always  aligns  with  the  specular  direction,  i.e.,  where 
6"  =  6‘ .  Conversely,  Figure  14  shows  that  for  the  2D  solution,  as  increases  the  SD 

peak  shifts  more  and  more  to  the  left  and  does  not  align  with  the  specular  direction^^. 

This  is  due  to  the  lack  of  the  aforementioned  amplitude  factor. 

Before  moving  on  in  the  analysis,  it  is  important  to  note  that  these  said  amplitude 
factors  are  also  responsible  for  the  discrepancy  seen  in  the  angular  SD  radius  (cf  the 
vertical  dashed  lines  in  Figure  10  and  Figure  13).  The  angular  SD  radius  only  comes 
close  to  the  1/e  locations  because  these  amplitude  factors  tend  to  “push”  out  the  wings  of 
the  far-field  power  distribution.  Nevertheless,  this  behavior  is  deterministic  in  nature  and 
the  angular  SD  radius,  as  derived  above,  adequately  characterizes  the  behavior  of  the 
scattered  SD. 


19.  The  full-wave  2D  MoM  solution  of  Basu  et  al.  shows  similar  behavior  [139,  140]. 
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5,2  Comparison  to  a  polarimetric  bidirectional  distribution  function 

The  analysis  presented  here  compares  the  3D  vector  solution  to  a  polarimetric 
bidirectional  distribution  function  (pBRDF)  developed  by  Priest  and  Meier  [49].  In 
general,  the  pBRDF  of  Priest  and  Meier  assumes  fully  incoherent  illumination;  thus,  there 
is  no  coherence  information  contained  within  the  solution.  Instead,  the  pBRDF  of  Priest 
and  Meier  provides  the  Mueller  matrix  for  statistically  rough  surfaces  that  are 
characterized  by  a  G-G  PDF  [cf.  Eq.  (138)];  as  a  result,  if  the  analysis  assumes 
unpolarized  illumination,  then  the  first  column  of  the  Mueller  matrix  (given  by  the 
pBRDF  solution)  becomes  the  scattered  Stokes  vector  [23,  30].  Provided  this  scattered 
Stokes  vector,  the  analysis  can  then  compare  the  normalized  scattered  SD  and  the 
scattered  DoP  between  the  pBRDF  solution  and  the  3D  vector  solution. 

5.2.1  Normalized  spectral  density  validation 

The  present  analysis  uses  the  Labsphere  Infragold  coupon  and  the  nominal  far- 
field  setup  (described  above).  It  also  assumes  unpolarized  illumination,  so  that 

and  =  0 .  Figure  15  shows  results  for  this  setup  with  partially  coherent 

illumination  at  non-normal  incidence,  where  =  0.5  and  O'  =  20° .  Note  the 

exact  agreement  between  the  pBRDF  and  3D  solutions  for  the  normalized  SD  with 
various  scattering  geometries.  Also  note  that  one  can  obtain  the  normalized  scattered  SD 
from  the  pBRDF  solution  by  cosine  correcting  the  first  term  of  the  scattered  Stokes 
vector  [34],  i.e.,  multiplying  by  cos  ) ,  and  dividing  by  the  max  value. 
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Figure  15.  Comparison  of  the  normalized  scattered  spectral  densities  obtained  from 
a  polarimetric  bidirectional  distribution  function  (pBRDF)  and  the  3D  vector 
solution  for  unpolarized  illumination  at  non-normal  incidence  and  a  very  rough 
conducting  surface,  (a)  depicts  an  in-plane  scattering  geometry,  whereas  (b)  depicts 
an  out-of-plane  scattering  geometry  with  results  as  a  function  of  a  single  polar  angle. 
Conversely,  (c)  and  (d)  depict  bi-static  scattering  geometries  as  a  function  of  a  single 
azimuth  angle.  Note  that  the  minimum  occurs  at  the  mono-static  observation  point 
in  both  (c)  and  (d). 
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5.2.2  Degree  of  polarization  validation 

The  initial  analysis  presented  here  again  uses  the  Labsphere  Infragold  eoupon  and 
the  nominal  far-field  setup  (deseribed  above).  It  also  assumes  unpolarized  illumination 
throughout,  so  that  and  =  0 .  Figure  16  shows  results  for  partially 

coherent  illumination  at  non-normal  incidence,  where  again  =  0.5  and  O'  =  20° . 

Note  the  exact  agreement  between  the  pBRDF  and  3D  solutions  for  the  scattered  DoP 
with  various  scattering  geometries. 

Now  the  analysis  assumes  a  tactical  engagement  scenario,  so  that  =1.064  pm , 
=  Tj  2  =  10  km ,  =  l.SAjl  cm ,  f  =  f  „„  =  0.25w^ ,  and  0'  =  56.4° .  This  setup 

corresponds  with  partially  coherent  illumination  of  NKB7  glass,  where  n  =  1.507  [143], 
at  Brewster’s  angle.  The  results  presented  in  Figure  17a  show  exact  agreement  between 
the  pBRDF  and  3D  solutions  for  scattered  DoP  assuming  very  rough  surface  conditions, 
where  I ^  =  lOO/ig  and  =  \0Aq  .  Figure  17b  then  shows  results  for  both  very  rough  and 
smooth-to-moderately  rough  surface  conditions  using  only  the  3D  vector  solution,  where 
again  f  ^  =  IOO/Iq  but  =  \0Aq  and  O.lA^ .  It  is  important  to  note  that  the  scattered  DoP 

only  exists  where  light  exists  in  the  analysis.  This  makes  sense  considering  that,  by 
definition,  the  scattered  DoP  depends  on  the  scattered  SD  [cf  Eqs.  (181)  and  (182)]. 

With  that  said.  Figure  17b  shows  that  the  scattered  DoP  does  not  depend  on  surface 
roughness,  at  least  for  the  assumptions  used  within  the  analysis  [cf  Footnote  10,  p.55].  It 
also  does  not  depend  on  coherence,  at  least  for  unpolarized  incident  illumination — this  is 
consistent  with  the  examples  given  by  Wolf  for  isotropic  beam  parameters  [13],  i.e., 
when  f  f  „„  (Appendix  E  shows  an  example  where  this  is  not  the  case). 
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Figure  16.  Comparison  of  the  scattered  degree  of  polarization  obtained  from  a 
polarimetric  bidirectional  distribution  function  (pBRDF)  and  the  3D  vector  solution 
for  unpolarized  illumination  at  non-normal  incidence  and  a  very  rough  conducting 
surface,  (a)  depicts  an  in-plane  scattering  geometry,  whereas  (b)  depicts  an  out  of 
plane  scattering  geometry  with  results  as  a  function  of  a  single  polar  angle. 
Conversely,  (c)  and  (d)  depict  bi-static  scattering  geometries  as  a  function  of  a  single 
azimuth  angle.  Note  that  the  minimum  occurs  at  the  mono-static  observation  point 
in  both  (c)  and  (d). 
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Figure  17.  Comparison  of  the  scattered  degree  of  polarization  for  unpolarized 
illumination  at  Brewster’s  angle  of  a  dielectric  surface  with  varying  roughness 
conditions,  (a)  shows  the  results  obtained  from  a  polarimetric  hidirectional 
distrihution  function  (pBRDF)  and  the  3D  vector  solution  for  non-normal  incidence 
at  Brewster’s  angle  with  very  rough  surface  conditions,  (h)  shows  the  results  from 
the  3D  vector  solution  for  both  very  rough  and  smooth-to-moderately  rough  surface 
conditions. 

5,3  Comparison  to  empirical  measurements 

To  compare  the  3D  vector  solution  with  empirical  measurements,  the  analysis 
uses  results  from  the  Complete  Angle  Scatter  Instrument  (CASI)  at  the  Air  Foree  Institute 
of  Teehnology  [146]  .  Figure  18  deseribes  the  seattering  geometry  assoeiated  with  the 
CASI  for  both  in-plane  measurements  (Figure  18a)  and  out-of-plane  measurements 
(Figure  18b).  Provided  Figure  18,  the  analysis  uses  the  following  angle  transformations, 
developed  by  Germer  and  Asmail  [148],  to  relate  the  CASI  scattering  geometry  to  that  of 
the  3D  veetor  solution  (of  Figure  9); 


O'  =  oos  ' 


oos(a")cos(y9*)J , 


(199) 


20.  Goldstein  readily  describes  the  principles  behind  dual-rotating-retarder  polarimety 
[147,  p.  357]. 
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(200) 


9"  =  cos 


-1 


cos  ( )  sin  (a* )  sin  ( )  +  cos  (a* )  cos  ( ji" )  cos  ) 


and 


'^TT 

f  =  tan-'  [b\a^)  +  —  -  tan^'  (/>' ,  a' ) . 


(201) 


In  Eq.  (201),  tan  ^  {b, a)  returns  the  inverse  tangent  of  b/a  after  taking  into  aeeount  the 


quadrant  of  the  point  {a,b) .  With  this  in  mind, 

¥  =  sin  )  j^eos  ( y' )  sin  {a' )  sin  )  +  cos  (sin  J 
+  cos  )  [^cos  (cr' )  cos  ( )  sin  )  -  sin  (cr' )  sin  ( y' )  J 

a  =-  cos  ( cos  ( y )  sin  ( ct' )  -  cos  ( a' )  sin  ( P' )  sin  ( y' ) 
+  sin ( S" ) [^cos [a' ) cos [y'  )- sin (ct' ) sin ( P' ) sin ( y' ) 

b‘  =  cos  {a' )  cos  ( y' )  sin  )  -  sin  (or' )  sin  ( r' )] , 

and  a  =-  cos  ( y" )  sin  (or' )  -  cos  (or' )  sin  (y^' )  sin  ( }^' ) . 


(202) 


(203) 


(204) 

(205) 


Provided  Eqs.  (199)-(205),  the  following  analysis  uses  the  Eabsphere  Infragold  eoupon 
and  the  nominal  far-field  setup  (described  above).  These  setup  parameters  best  mateh 
those  used  by  the  CASE  It  also  assumes  unpolarized  partially  coherent  illumination  at 
non-normal  incidence,  so  that  A^  =  A^,  B^=0 ,  or^  =  or^„  =  0.5 ,  and  O'  =  20° . 

5.3.1  In-plane  measurements 

Eor  in-plane  measurements,  P'  =  y'  =  0  .  Asa  result,  the  angle  transformations 

given  above  simplify,  so  that  or'  =9'  =  20°  and  S"  =  9" .  Eigure  19  shows  results  for  the 
in-plane  comparison  study.  Note  that  the  results  for  the  CASI  measurements  are  highly 
oseillatory  for  both  the  normalized  seattered  SD  and  the  scattered  DoP.  This  is  most 
likely  due  to  the  faet  the  measurements  eontain  speckle.  Spinning  the  sample  would 
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average  out  this  speekle;  however,  the  eurrent  CASI  setup  does  not  allow  for  eontinuous 
spinning  of  the  sample.  In  addition,  the  CASI  measurements,  with  respeet  to  the 
normalized  seattered  SD,  show  that  the  Labsphere  Infragold  eoupon  does  not  maintain  a 
G-G  model  for  the  underlying  surfaee  statisties.  Instead,  a  stretehed  exponential- 
stretehed  exponential  (SE-SE)  model  better  eharaeterizes  the  Eabsphere  Infragold  eoupon 
[138].  This  is  seen  by  eomparing  the  results  to  a  full-wave  2D  MoM  solution  [139,  140]. 
Here,  the  analysis  obtained  the  seattered  field  from  400  independent  rough  surface 
realizations  simulated  using  the  method  described  by  Yura  and  Hanson  [145]  with  both 
G-G  and  SE-SE  PDEs  [cf.  Eq.  (138)].  Before  moving  on  in  the  analysis,  it  is  important  to 
note  that  the  CASI  measurements  do  not  contain  enough  fidelity,  with  respect  to  the 
scattered  DoP,  to  compare  them  to  those  obtained  by  the  3D  vector  solution. 

5.3.2  Out-of-plane  measurements 

Eor  the  out-of-plane  measurements,  the  CASI  used  the  following  setup 
parameters:  a'  =  20°  and  =  25° .  Accordingly,  one  can  use  the  angle  transformations 
given  above  in  Eqs.  (199)-(205)  to  relate  the  3D  vector  solution  to  the  CASI 
measurements.  Eigure  20  shows  results  for  the  out-of-plane  comparison  study.  Note 
again  that  the  CASI  measurements,  with  respect  to  the  normalized  scattered  SD,  show 
that  the  Eabsphere  Infragold  coupon  does  not  maintain  a  G-G  model  for  the  underlying 
surface  statistics.  Also  note  that  the  CASI  measurements  do  not  contain  enough  fidelity, 
with  respect  to  the  scattered  DoP,  to  compare  them  to  those  obtained  by  the  3D  vector 
solution. 
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(a) 


Detector 


(b) 


Figure  18.  Description  of  the  Complete  Angle  Scatter  Instrument  at  the  Air  Force 
Institute  of  Technology,  (a)  shows  the  scattering  geometry  used  to  collect  in-plane 
measurements  [140],  whereas  (h)  shows  the  scattering  geometry  used  to  collect  out- 
of-plane  measurements  [148]. 
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Figure  19.  Comparison  between  in-plane  measurements  obtained  with  tbe 
Complete  Angle  Scatter  Instrument  (CASI)  at  tbe  Air  force  Institute  of  Tecbnology, 
a  full  -wave  2D  method  of  moments  (MoM)  solution  with  two  different  surface 
models,  and  tbe  3D  vector  solution  for  unpolarized  illumination  at  non-normal 
incidence  of  a  very  rough  conducting  surface,  (a)  and  (c)  show  the  normalized 
scattered  spectral  density  and  the  scattered  degree  of  polarization  as  a  function  of  a 
single  polar  angle,  (b)  and  (d)  show  the  same  results  in  log  scale. 
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Figure  20.  Comparison  between  out-of-plane  measurements  obtained  with  tbe 
Complete  Angle  Scatter  Instrument  (CASI)  at  tbe  Air  force  Institute  of  Tecbnology 
and  tbe  3D  vector  solution  for  unpolarized  illumination  at  non-normal  incidence  of 
a  very  rough  conducting  surface,  (a)  and  (c)  show  tbe  normalized  scattered  spectral 
density  and  tbe  scattered  degree  of  polarization  as  a  function  of  a  single 
transformation  polar  angle,  (b)  and  (d)  show  tbe  same  results  in  log  scale. 
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5.4  Comparison  to  a  paraxial  solution 

An  individual  can  analyze  the  3D  vector  solution  within  a  Cartesian  coordinate 
system.  With  that  said,  one  must  use  the  following  relationships: 


Provided  the  relationships  in  Eq.  (206),  the  analysis  can  compare  the  3D  vector  solution 
to  an  ABCD-matrix  approach  outlined  by  Korotkova  [115].  In  general,  the  ABCD- 
matrix  approach  describes  paraxial  wave  propagation  through  any  complex  optieal 
system.  When  modeling  rough  surface  scattering  using  the  ABCD-matrix  approach,  a 
phase-screen  transmittance  function  characterizes  the  scattering  proeess  for  very  rough 
surface  conditions.  The  inclusion  of  a  soft-Gaussian  aperture  in  the  model  also  accounts 
for  the  size  of  the  scattering  surface  and  accompanying  diffraction  effects. 

The  present  analysis  assumes  horizontally  polarized  (s-pol)  illumination  and  a 
taetical  engagement  scenario,  so  that  =  0  ,  /if,  =  1 .064  pm , 


:  10  km,  =  2.54/2  cm. 


:  0.25w  ,  and  0' 


This  setup 


corresponds  with  partially  coherent  illumination  at  normal  incidence.  To  ensure  that 
incident  illumination  “fits”  on  the  rough  surface,  the  present  analysis  removed  the  soft- 
Gaussian  aperture  in  the  ABCD  solution.  In  addition,  the  present  analysis  assumed  an 
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idealized  perfeetly  refleeting  rough  surfaee,  where  =  lOO/lo  and  =  IO/Iq  .  Sueh 

surfaee  statisties  correspond  to  very  rough  surface  conditions  [cf  Eq.  (174)]. 

Provided  this  setup,  Figure  21  shows  excellent  agreement  between  the  ABCD 
paraxial  solution  and  the  3D  vector  solution.  For  instance,  the  magnitude  of  the  scattered 
SDoC  (Figure  21a)  shows  exact  agreement.  The  vertical  line  depicted  here  is  the  angular 
SDoC  multiplied  by  the  propagation  distance  to  observation.  Note  that  it  correctly 
identifles  the  1/e  location.  In  addition,  the  normalized  scattered  SD  shows  relatively 
good  agreement  between  the  two  solutions.  The  wings  of  the  power  distribution 
associated  with  the  3D  vector  solution  tend  to  be  wider  than  those  associated  with  the 
ABCD  solution;  nonetheless,  the  results  show  good  agreement  out  to  the  vertical  lines 
depicted  in  Figure  21b.  These  vertical  lines  are  the  angular  SD  radius  multiplied  by  the 
propagation  distance  to  observation.  Furthermore,  the  scattered  DoP  (Figure  21c)  shows 
exact  agreement  between  the  two  solutions.  Based  on  the  assumptions  used  within  the 
analysis  [cf.  Footnote  10,  p.55],  no  de -polarization  occurs  upon  scattering  from  the 
perfectly  reflecting  rough  surface  [22],  at  least  for  isotropic  beam  parameters,  i.e.,  when 
^ XX -^uu  (Appendix  E  shows  an  example  where  this  is  not  the  case).  This  is  also  the  case 

for  p-pol  illumination  at  normal  incidence  and  is  consistent  with  the  2D  scalar-equivalent 
solution  [10,  11]. 
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Figure  21.  Comparison  between  an  ABCD  paraxial  solution  and  the  3D  vector 
solution  for  partially  coherent  illumination  at  normal  incidence  of  a  very  rough 
perfectly  reflecting  surface,  (a)  shows  the  magnitude  of  the  scattered  spectral 
degree  of  coherence  as  a  function  of  the  distance  between  two  values  in  the  x 
direction,  whereas  (b)  shows  the  scattered  normalized  spectral  density  and  (c)  shows 
the  scattered  degree  of  polarization  as  a  function  of  a  single  value  in  the  x  direction. 
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6  Conclusion 


This  dissertation  develops  a  3D  veetor  solution  for  the  far-field  seattering  of 
spatially  partially  eoherent  eleetromagnetie  beam  illumination  from  statistieally  rough 
surfaees.  Compared  to  previous  researeh  efforts,  the  3D  veetor  solution  developed  in  this 
dissertation  signifieantly  extends  the  rough  surfaee  seattering  literature,  sinee  most  of  the 
theory  developed  to  date  eonsiders  only  the  effeets  of  fully  eoherent  and  fully  ineoherent 
illumination  in  the  formation  of  solutions.  Moreover,  the  analysis  contained  in  this 
dissertation  considers  three  different  material  substrates:  dielectrics,  conductors,  and  a 
perfect  electrical  conductor  (PEC).  By  incorporating  the  effects  of  the  material 
parameters  within  the  analysis,  the  3D  vector  solution  developed  here  transcends 
previous  efforts  which  included  the  effects  of  partially  coherent  beam  illumination,  but 
not  the  effects  of  the  material  parameters. 

To  develop  the  3D  vector  solution  contained  in  this  dissertation,  the  analysis  uses 
the  physical  optics  approximation  (Kirchhoff  boundary  conditions).  It  also  uses  a 
Gaussian  Schell-Model  (GSM)  form  for  the  incident  field  cross  spectral  density  matrix 
(CSDM).  This  allows  for  the  formulation  of  a  closed-form  expression  for  the  scattered 
field  CSDM.  In  practice,  two  separate  analytical  forms  result  for  the  scattered  field 
CSDM.  The  first  analytical  form  is  applicable  to  smooth-to-moderately  rough  surfaces. 

It  depends  on  an  infinite  series  that  is  slowly  convergent;  however,  one  can  still  gain 
physical  intuition  from  the  resultant  closed-form  expression.  The  second  analytical  form 
is  applicable  to  very  rough  surfaces.  Here,  the  closed-form  expression  is  extremely 
physical,  and  under  certain  circumstances,  maintains  a  GSM  form.  This  is  in  agreement 
with  published  results  valid  only  in  the  paraxial  regime. 
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By  using  the  CSDM  notation  within  the  3D  veetor  solution,  all  aspeets  inherent  in 
Wolfs  unified  theory  of  eoherenee  and  polarization  apply  here  [12,  13].  Aeeordingly, 
the  analysis  eontained  in  this  dissertation  readily  formulates  the  speetral  degree  of 
eoherenee  (SDoC),  the  speetral  density  (SD),  and  degree  of  polarization  (DoP)  to  observe 
the  eoherenee,  size,  and  polarization  properties  of  the  seattered  radiation  in  eomparison  to 
the  ineident  radiation.  The  SDoC,  SD,  and  DoP  are  measurable  quantities  in  praetiee  and 
serve  as  metrics  in  which  to  compare  the  3D  vector  solution  to  previously  validated 
solutions  and  empirical  measurements.  With  this  said,  the  analysis  contained  in  this 
dissertation  rigorously  compares  the  obtained  3D  vector  solution  to  previously  validated 
solutions  and  empirical  measurements.  The  comparisons  show  good  agreement  under  the 
appropriate  conditions. 

6,1  Contributions 

There  are  two  significant  contributions  which  have  already  resulted  from  this 
dissertation  effort  (with  more  to  come).  The  first  is  a  SPIE  Newsroom  article  which 
highlights  the  contributions  of  the  2D  scalar-equivalent  solution  [10,  11]  and  sets  the 
stage  for  the  3D  vector  solution  developed  here  [8].  This  article  also  proposes  a  novel 
experiment,  which  could  validate  the  analysis  contained  in  both  the  2D  and  3D  solutions 
and  is  an  area  for  future  research  (discussed  below).  The  second  is  a  16*  Annual 
Directed  Energy  Symposium  proceeding  which  formulates  the  3D  vector  solution  for  a 
PEC  material  substrate  [9].  This  contribution  was  selected  for  publication  in  both  the 
2013  AP-URSI  and  2014  IEEE -Aerospace  conferences  based  on  abstract  submission; 
however,  due  to  setbacks  in  the  federal  government’s  budget,  this  publication  was 
postponed. 
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6,2  Future  areas  of  research 


Future  work  in  rough  surface  scattering  using  partially  coherent  illumination 
includes,  but  is  not  limited  to,  the  following  topics. 

•  Anisotropic  material  substrates,  where  the  surface  correlation  length  is 
directionally  dependent  within  the  analysis  [43,  66]. 

•  Partially  coherent  pulsed  illumination  that  emanates  from  a  source  that  is  not 
statistically  wide-sense  stationary  [82,  149]. 

•  The  construction  of  designer  surfaces  which  control  the  spatial  and  temporal 
coherence  properties  of  the  scattered  radiation  [85,  113]. 

•  The  use  of  a  non  uniformily  polarized  source,  where  the  source  width  is 
directional  dependent  within  the  analysis  [101,  135]. 

•  Surfaces  with  large  surface  slopes  so  that  shadowing,  masking,  and  multiple 
scattering  effects  play  a  role  in  the  analysis  [66,  150]. 

•  The  development  of  a  closed-form  coherent  polarimetric  bidirectional  distribution 
function,  which  is  validated  using  high-fidelity  measurements  [46,  127]. 

•  The  development  of  wave-optics  simulations  to  model  the  spatially  partially 
coherent  beam  illumination  and  the  laser-target  interaction  [102,  103]. 

•  Validation  experiments  which  use  a  broadband  source  and  narrowband  filters  to 
look  at  changes  in  the  SD  and  DoP  upon  scattering  off  a  rough  surface  [8,  13]. 
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Appendix  A.  Using  the  geometrical  optics  approximation  to  relate  the 

tangential  fields  at  the  scattering  surface 

Referencing  the  micro-scale  scattering  geometry  given  in  Figure  23  [18],  the  law 
of  reflection  dictates  that  the  angle  of  incidence  equals  the  angle  of  reflection,  . 

This  is  a  direct  result  from  geometrical-optics  (GO)  approximation  [22],  As  such,  the 
following  relationships  hold  true  at  the  closed  surface  S'  of  a  scattering  object: 


(207) 

h-k“'^=-cos(r^). 

(208) 

h-k^‘^^=cos(^^'^^). 

(209) 

hxk'""=hxk''‘'-^. 

(210) 

hxk'^=tsin(^'""). 

(211) 

hxk^‘'^=tsin(^''"^). 

(212) 

where  h  is  the  unit  normal  vector,  k'”''  is  the  unit  incident  propagation  vector,  is  the 
unit  reflected  propagation  vector,  and  t  is  the  unit  tangential  vector.  Provided  Eqs.  (207) 
-(212),  the  purpose  of  this  appendix  is  to  demonstrate  how  the  tangential  reflected  fields, 
hxE''”'^  and  relate  to  the  tangential  incident  fields,  hxE'""  and  hxH'“,  for 

different  material  substrates,  i.e.,  dielectrics,  conductors,  and  a  perfect  electrical 
conductor  (PEC). 
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Substrate 


Figure  22:  The  micro-scale  geometry  describing  how  to  relate  the  incident  field  to 
the  refiected  field  using  the  geometrical  optics  approximation. 

The  following  relationship  defines  the  Fresnel  reflection  coefficient  at  the 
closed  surface  S\ 


rs  = 


E" 


E 


ref 


(213) 


so  that  h  X  =  r^h  x  E'"" .  (214) 

Equation  (31)  above  repeats  this  GO  result.  Moreover,  since  [14] 

E  =  -;7okxH,  (215) 

a  similar  relationship  exists  for  the  magnetic  fields  to  that  given  in  Eq.  (214)  for  the 
electric  fields.  Substituting  Eq.  (215)  into  Eq.  (214), 

-//oh  X  X  )  =  -//or^h X  (k'""  X  H'“ ) .  (216) 

Eurthermore,  using  the  following  vector  identity: 

Ax(BxC)  =  B(A-C)-C(A-B)  (217) 

and  Eqs.  (207)-(212),  Eq.  (216)  manipulates  as  follows: 
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^  (n  ■  (n  ■  (n  ■  H'"" )  -  r^H'""  (n  ■  k“" ) 

^  (n  ■  H'"" )  +  (n  ■  k'"'^ )  =  r^k'”"  (n  ■  H'"" )  -  r^H'^  (n  ■  k“" ) 

^  r> X  k'-^^  (n  ■  H'"" )  +  n  X  (n  ■  k“" )  =  r> x  k'""  (n  ■  H'"'^ )  -  r^n x  H'""  (n  ■  k“'^ ) 

^  r> X  k'""  (n  ■  )  +  n X  (n  ■  k“" )  =  r^n  x  k'"'^  (n  ■  )  -  r^n  x  H'"'^  (n  ■  k“" ) 


.(218) 


From  Eq.  (218),  it  is  readily  apparent  that 

hxH'■‘^^=-r5hxH"'^  (219) 

Equation  (32)  above  repeats  this  GO  result. 

The  tangential  eomponents  of  the  eleetrie  field  equate  zero  along  the  surface  of  a 
PEC  [1].  Thus,  the  following  relationship  holds  true: 

h  X  (e'""  +  )  =  0  ^  h  X  E'"'^  =  -h  X  .  (220) 


Equation  (37)  above  repeats  this  GO  result.  Substituting  Eq.  (215)  into  Eq.  (220), 

//oh  X  (k'”"  X  H'“ )  =  -//oh  X  X  ) .  (221) 

Additionally,  using  Eq.  (217)  and  Eqs.  (207)-(212),  Eq.  (221)  manipulates  as  follows: 

^  k-  (n .  )  -  H'”"  (h  ■  k“" )  =  -  (h  ■  )  -  H''^^  (h  ■  k"*^^ ) 

^  k'-  (h  ■  H''^^ )  -  H'""  (h  ■  k'''^-^ )  =  (h  ■  H''"^ )  -  (h  ■  ) 

^  h  X  k'"'^  (h  ■  )  -  h  X  H'""  (h  ■  k""^^ )  =  h  X  k''^^  (h  ■  )  -  h  x  (h  ■  k'''^^ ) 

^  h X  k'-^^  (h  ■  )  -  h  X  (h  ■  k"^^ )  =  h  X k'-"^  (h  ■  )  -  h  X  (h  ■  k'''^-^ ) 


.  (222) 


Erom  Eq.  (222),  it  is  readily  apparent  that 

hxH''‘'^=hxH“C  (223) 

Equation  (38)  above  repeats  this  GO  result. 
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Appendix  B.  Using  the  physical  optics  approximation  to  simplify  the 


equivalent  surface  current  densities 


The  physical-optics  (PO)  approximation  uses  the  geometrical-optics  (GO) 
approximation  to  formulate  the  current  densities  involved  in  an  electromagnetic 
scattering  problem.  With  that  said,  the  purpose  of  this  appendix  is  to  demonstrate  how 
the  equivalent  surface  current  densities,  which  result  from  using  surface  equivalence, 
simplify  to  known  quantities  when  using  the  GO  approximation  in  association  with  the 
PO  approximation. 

To  make  this  concept  manifest,  first  consider  the  micro-scale  scattering  geometry 
given  in  Figure  23.  Here,  the  law  of  reflection  dictates  that  the  angle  of  incidence  equals 
the  angle  of  reflection,  .  This  is  a  direct  result  from  the  GO  approximation 

[22].  In  using  the  GO  approximation  as  part  of  the  PO  approximation,  one  must  assume 
that  the  scattering  object  and  its  associated  surface  curvature  are  large  compared  to  the 
wavelength  of  the  incident  fields.  Such  approximations  are  exact  when  the  scattering 
surface  is  homogeneous,  infinite,  and  planar  [14,  22].  As  such,  when  using  the  PO 
approximation,  the  following  relationships  hold  true  at  the  closed  surface  5*  of  a  scattering 
object: 


s  =  -^ - -  s  =  k“'^xp  =  ,  (224) 

k'""xn 

p  =  sxk“"  -p  =  k“"xs,  (225) 

=  -s  X  k'-^^  X  s ,  (226) 

and  -h-k“"=h-k''"^,  (227) 
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where  s  is  the  unit  perpendicular  vector,  p  is  the  unit  parallel  vector,  is  the  unit 
parallel  vector  upon  reflection,  h  is  the  unit  normal  vector,  k'"'^  is  the  unit  incident 
propagation  vector,  and  is  the  unit  reflected  propagation  vector.  Provided  Eqs.  (224) 
-(227),  one  can  then  write  the  equivalent  surface  current  densities,  and  ,  in  terms 

of  only  the  known  incident  electric  field  E'""^ .  This  is  shown  in  the  analysis  to  follow. 


Figure  23:  The  micro-scale  geometry  describing  how  to  relate  the  incident  field  to 
the  refiected  and  transmitted  fields  using  the  geometrical  optics  approximation. 

As  stated  in  Section  2.3,  when  using  the  PO  approximation  with  surface 
equivalence,  one  replaces  the  unknown  scattered  fields  with  reflected  fields.  In 
particular,  ~  and  H*"'  ~  ,  so  that  the  equivalent  surface  current  densities, 

J|'  and  become 

=  h  X  =  h  X  ( H'""  +  H""' )  -  h  X  ( )  (228) 

and  m;''  =  -h  X  E“'  =  -h  x  (e''“  +  E""'  )  -  -h  x  (e“"  +  E'''^^  ) .  (229) 

Moreover,  in  using  the  following  relationship  [14]; 
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(230) 


H  =  — kxE 


and  writing  the  respective  electric  fields  in  terms  of  their  perpendicular  and  parallel 
components,  i.e., 


E  =  (ss  +  pp)-E , 

(231) 

the  equivalent  surface  current  densities,  and  ,  become 

^ h X  Tk'""  X  (ss  +  pp )  ■  E“"  +  k''"-^  (ss  +  ’ 

Vo 

(232) 

and 

=-hx[(ss  +  pp)-E'"'^  +  (ss  +  p^‘'^p''‘^-^)-E''‘^-^" . 

(233) 

Defining  the  perpendicular  and  parallel  Fresnel  reflection  coefficients  as 

s-E''*'^  ;7cos(^“")-;7ocos(^"'") 

'2  -  g  .  ^inc  -  ^  J  _  J 

(234) 

and 

^ref  ^ref  COS  (6>'^  )  +  7/  COS  [O"''  ) 

“  p  ■  E'^  “  cos  ( )  + ;/  cos  ( )  ’ 

(235) 

where  is  the  angle  of  transmission,  is  the  impedance  of  free-space,  and  rj  is  the 
impedance  of  the  material  substrate  (cf  Figure  23),  one  can  then  write  the  equivalent 
surface  current  densities,  and  ,  in  terms  of  only  the  known  incident  electric  field 

E'“.  Specifically, 

=  — hx[k'"^x(ss  +  pp)  +  k^‘'^x(srJ  +  p^‘^^ijip)]-E“"  (236) 

Vo 

and  M5^=-hx  (ss  +  pp)  +  (sr^s  +  p''‘''^i]ip)J-E“‘' .  (237) 
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Before  moving  on  in  the  analysis,  it  is  important  to  note  that  the  Fresnel  refleetion 
eoeffieients,  as  defined  in  Eqs.  (234)  and  (235),  are  readily  derived  throughout  the 
eleetromagnetic  and  op  ties  literature.  Balanis’  treatment  is  partieularly  insightful  [14],  as 
his  problem  setup  direetly  relates  to  that  used  in  Figure  23. 

Using  the  relationships  found  in  Eqs.  (224)-(227),  the  equivalent  surfaee  current 
densities,  and  as  given  in  Eqs.  (236)  and  (237),  simplify  even  further,  such  that 


and 


r/=_hx 

Vo 


X  (ss  +  pp)  +  X  (sr^s  +  )1  ■  E” 


=  —  h  X  fsp  -  ps  +  p''‘''^r|  s  -  si:p]  •  E'"" 

Vo  ^  ^ 

lnx[s(l-^,)p-(p-p"^rqs].E'"' 

hxs(l-ij|)p-(hxsxk'”‘'  -hx-sxk''‘''^rj^  js 


Vo 
1 

Vo^ 


■E” 


M^'*'  =  -h  X  j^ss  +  pp  +  srj^s  +  p''‘'^i]|p]  ■  E'“ 

]s(l  +  rJs  +  (p  +  r"'i,)p]^E'“ 

nxs(l  +  r^)s+(nxp  +  nxp'^^/ji)p1-E'"' 


=  -nx 


Eurthermore,  using  the  following  vector  identity: 

Ax(BxC)  =  B(A-C)-C(A-B) 


(238) 


(239) 


(240) 


one  can  then  write  the  equivalent  surface  current  densities,  and  ,  in  terms  of  the 

dyadics,  and  ,  where 
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and 


jeq 

S 


1  - 

1 

^0 


nxs(l-?]|)p-(nxsxk'"''  -nx-sxk''‘'^r^  js 


(« 

^inc 

)+S(n 

^ref 

i  1 

(« 

^inc 

)-4" 

^inc 

) 

1 - 1 

s )  ■  E'”'" ,  since  n  ■  s  =  0 

g,  g,«c 


nxs(l-?]|)p-s(n  ■k'"‘')(l-r^)s 
— nxj^s(l-?]|)p-sxk“‘'  (l-ri)s 


E" 


(241) 


=  — nx[s(l-?]|)p-p(l-r^)s 
^0 


■E” 


^0 


m;^  =  - 


=  -nx 


nxs(l  +  r^)s  +  (nxp  +  nxp''‘'^?]i)pJ-E“‘' 
nxs(l  +  r^)s  +  (nxsxk'"''  +nx-sxk''‘''^?]i)p 
nxs(l  +  r^)s+  s(n  -k'""  j-s(n  ■k''*"'^  j?]iP 
nxs(l  +  r^)s+  s(n-  k'"'" )  +  s  (n  ■  k'""  j  ?]|P 
nxs(l  +  r^)s  +  s(n-  k'"^'  j  (l  +  /]|  )p 
s(l  +  r^)s  +  sxk'“  (l  +  ?]|)pJ-E 
s(l  +  rJs  +  p(H-,i)p]-E"' 


■  E'”'" ,  since  n  ■  s  =  0 

l^mc 


(242) 


■  E" 


:-nx 


_M^?.E'- 


Equations  (1 15)-(1 18)  above  repeat  these  PO  results. 


no 


Appendix  C.  Using  the  method  of  stationary  phase  to  solve  the  integrals 
with  respect  to  the  plane-wave  spectrum  representation 

The  plane-wave  speetrum  representation  often  results  in  rather  eomplex  integral 
expressions.  In  praetiee,  one  must  employ  mathematieal  teehniques,  sueh  as  the  method 
of  stationary  phase  (MoSP),  to  solve  these  eomplex  integral  expressions.  The  purpose  of 
this  appendix  is  to  demonstrate  how  one  uses  the  MoSP  to  reduee  the  integrals  with 
respect  the  plane-wave  spectrum  representation  into  closed-form  expressions. 

Within  the  analysis  given  above,  the  plane-wave  spectrum  representation  results 
in  integral  expressions  which  take  the  following  form  [cf.  Eqs.  (141)-(144)]; 

(243) 

where  f  =  f  ^2 ■> Tp T2 )  is  a  representative  amplitude  dyadic  that 

contains  all  of  the  amplitude  terms  and  g  =  g  ( ’  -^2  ’  Ti  ’  T2 )  is  ^  phase 

function  that  contains  all  of  the  phase  terms.  In  using  the  MoSP  (cf  Section  2.4.3),  Eq. 
(243)  simplifies  into  the  following  closed-form  expression; 

^2;7-V  f  f  TT  \ 

- r  °  ,|i/2  ^xp[Ago]gxp  7-sgn{a"go}  (244) 

Uoj  |Det{a^g4|  V  4 

where  fg  and  gg  are  the  amplitude  dyadic  and  phase  function  evaluated  at  the  critical 
points  of  the  first  kind,  respectively.  In  general. 
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9"g 

d^g 

d^g 

k:^K2 

KxK2 

y-g 

a^g 

^K2^Kx 

^K2^K2 

yg 

a^g 

KrK. 

KM2 

KxK2 

9"g 

9"g 

y-g 

9"g 

K2K. 

K2^K2 

K2KX 

HKif 

so  that  d^gg  is  a  real  symmetric  non-degenerate  matrix  which  is  also  evaluated  at  the 

critical  points  of  the  first  kind.  To  determine  the  critical  points  of  the  first  kind,  one  must 
first  account  for  the  phase  function  g  . 

In  the  analysis  given  above,  the  following  common  phase  function  results  [cf.  Eq. 


(149)]; 


g  =  [r,+{y-y)y2~\K^-[r^+{y-y)y[~\K 


-I- — 


^  [AKi  - - y'lKi ) : 


1  \  A  XA  1  XI/  j 

Kq  Kq 


+  (x-fX-x-f2x;)  +  (y-fjy;-y-f2y;) 


(246) 


where 


Ki2  = 


1- 


^  k‘  ^ 


^  rJ  \ 
^ul,2 


V  -0  y 


V  -0  y 


(247) 


Provided  Eq.  (246),  one  then  computes  the  following  first-order  partial  derivatives: 


dg  _  x;  ,  Kx[r,+{y-x)y[~\ 


dkt 


1- 


V'"o  y 


f  ri  y 

^ul 


V  '"0  y 


(248) 
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dg  _X2 

dk[,  k. 


Kil^s  +  iy  ■^)y2~\ 


1 


f  P  ^ 

^x2 

\K  j 


\K  j 


(249) 


(y-»)X  ^  Kx\r,+{y-x)y[~\ 


1 


1- 


^x\ 

\K  J 


f  ij  \ 


\K ) 


(250) 


and 


dg  _(y'u)j;  ^L[^.+(y-v)j2] 


dk. 


u2 


\K  J 


f  \ 

^u2 

\K  J 


(251) 


By  setting  Eqs.  (248)-(251)  equal  to  zero  and  solving  for  A:(j  j  and  k^^  j ,  one  then  obtains 


four  equations  and  four  unknowns.  Speeifieally, 


9g 


dk 


=  0  ^  Ki,i  =  ± 


y,2'\j^0  {j^u\,2 ) 


xl,2 


+  <2  +  (y  ■  v)  y'2  +  2r,  (y  ■  v)  X2 


(252) 


and 


dg 


dk. 


=  0: 


k’  -  + 
%1,2  — 


(y-«)yi' 


k 

2\l'^0 


E;.4 


uL2 


+  (y  ■ « )"  Eu  +  (y  •  v)'  X 2  +  (y  ■  v)  X: 


2  /2 


(253) 


These  four  unknowns  are  the  eritieal  points  of  the  first  kind,  i.e.,  Xi  2  and  k[^  2  •  ^y 

substituting  Eq.  (252)  into  Eq.  (253)  and  viee  versa,  one  then  determines  the  eritieal 
points  of  the  first  kind  as 


J^l  _  '  1^1 

^x\.2  *^1,2  ^u\,2 


^o(y-u)  / 

- Tu  • 


(254) 


Equation  (156)  above  repeats  this  MoSP  result. 
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In  writing  Eq.  (254),  one  replaces  ±|x'2|  and  ±|ei2|  (which  result  from  solving 
quadratic  expressions)  with  x'2  and  2 .  This  result  is  physically  intuitive  based  on  the 
macro-scale  scattering  geometry  given  in  Figure  9a.  Furthermore,  one  assumes  that 

(255) 

This  is  analogous  to  saying  that  distance  from  source-plane  origin  to  the  surface-plane 
origin  is  much  greater  than  half  the  surface  length,  i.e.,  »  F  ,  which  is  typically  the 

case  for  most  directed-energy  and  remote-sensing  engagement  scenarios.  To  provide 
some  idea  of  how  much  greater,  letting  r^=5L  ,  r^=  25 L  ,  and  =  lOOF  results  in 
percentage  errors  of  18%  ,  4% ,  and  1%  ,  respectively.  These  percentage  errors  result 
from  setting  Xj2  =7;  2  =  L,  y  •  v  =  sin(^' )  =  1 ,  and  y  •  v  = -cos(^' )  =  0  inFq.  (255). 

Provided  Fqs.  (248)-(251),  one  then  computes  the  following  second-order  partial 
derivatives: 


(^xif  [F  +  (y-v)yi] 


.+(y-v)y; 


(256) 


(^;ii)'[F+(y-v)X]  ^ 


F+(y-v)Ti 


(258) 
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[K2)\r,+{y-y)y[~\ 


';+(y-v)j; 


^[K2) 


fij  A 


x2 


k‘ 


kK  j 


f  i^i  A 

^Ul 

kK  j 


-|3/2 


kiw- 


f  U  \ 

^x2 

kK  j 


^  t'2  ^ 
kK  j 


(259) 


5^g 


5^g 


dk'M,  dki^ki,  dki^ki,  dk:M2 


=  0, 


(260) 


d^g  _  KK\rs+{y-'^)y['\ 


^xl 


\^0  J 


f  li  ^ 

^u\ 


\^0  J 


3/2 


0, 


(261) 


and 


5^g 


K2K2[r,  +  {y-y)y2~\ 


K2K2 


V2 

\K  j 


\K  j 


3/2 


(262) 


In  Eqs.  (261)  and  (262),  one  assumes  that  k^  »  A:/  and  k^  »  A:^  (reasons  diseussed 
below);  eonsequently,  d^g  beeomes  a  diagonal  matrix  upon  substituting  Eqs.  (256)-(262) 
into  Eq.  (245).  Subsequent  evaluation  of  d^g  at  the  eritical  points  of  the  first  kind  [cf. 

Eq.  (254)]  results  in  the  following  relationships: 


Det{a"go}  =  - 


[fi  +  (y  •  v)  El']"  [r,  +  (y  •  v)  y'J 


^r+(y-v)  yr-r; 


^"  +  (y-v)  E--r/ 


(263) 


and  sgn{a^go}  =  /l^{a^go}-2  {a^go}-0,  (264) 

where  2+  {3^go}  are  the  number  of  positive  and  negative  eigenvalues  of  d^gg .  Equations 
(152)-(155)  above  use  these  relationships  in  eombination  with  Eq.  (244). 
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As  discussed  above  in  Seetion  4.3,  the  analysis  assumes  that  »  k[  and 
kl'»  kl-  In  so  doing,  the  analysis  also  assumes  that 


k 


0’ 


inf 


ing 


(265) 


Thus,  in  evaluating  the  phase  funetion  g  [ef.  Eq.  (246)]  at  the  eritieal  points  of  the  first 
kind  [ef  Eq.  (254)],  the  following  relationship  results: 


go 


X. 


n 


(y-u) 


2  /2 

Ti 


2r  2r 


2r 


2r 


+(y-v)y2-(y-v)yi' 


(266) 


+  (x-fX-x-f2x;)  +  (y-fiy;-y-f2y;) 


Equations  (152)-(155)  and  subsequently  Eq.  (157)  above  use  this  relationship  in 
eombination  with  Eq.  (244).  With  this  in  mind,  the  analysis  evaluates  the  representative 

amplitude  dyadie  f  [ef  Eqs.  (145)-(148)]  at  the  eritieal  points  of  the  first  kind  below  in 
Appendix  D. 
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Appendix  D.  Defining  the  dyadics  that  contain  all  of  the  amplitude  and 

phase  terms  evaluated  at  the  critical  points  of  the  first  kind 

In  using  the  MoSP  to  simplify  Eqs.  (141)-(144),  Eqs.  (152)-(156)  result.  Inherent 
in  Eq.  (152)-(155)  are  the  following  dyadies: 


L  L  L  L 


jC  =  1  1  1  1  f/  dx[dx2dy[dy2 , 

-L  -L  -L  -L 

(267) 

L  L  L  L 

/c  =  j  j  j  j  dx[dx^dy[dy[  , 

-L-L  -L  -L 

(268) 

L  L  L  L 

jCi  =1111  dx\dxjdy\dy^  , 

-L  -L  -L  -L 

(269) 

and 

L  L  L  L 

AT  =  1  1  1  1  dx[dx2dy[dy2 , 

-L  -L  -L  -L 

(270) 

where  if ,  if ,  ,  and  are  the  amplitude  dyadics  evaluated  at  the  critical  points  of 

the  first  kind,  respectively,  and  gg  is  a  common  phase  function  which  is  also  evaluated  at 
the  critical  points  of  the  first  kind.  Provided  Eqs.  (267)-(270),  the  purpose  of  this 

appendix  is  to  explicitly  define  the  dyadics  C,  K,  Ad ,  and  Af  for  different  material 
substrates,  i.e.,  dielectrics,  conductors,  and  a  perfect  electrical  conductor  (PEC). 

General  problem  for  dielectrics 

With  some  mathematical  prowess,  one  establishes  the  dyadics  given  in  Eqs.  (267) 
-(270)  as 
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and 


'^xx 

r 

'^xy 

r 

'^xz 

4. 

4 

4. 

=  4j  (i  =  x,y,z;  j  =  x,y,z} 


a:= 


AjAj 

m=x. 

1  ^ 

,w  n=x,j 

4. 

4. 

4. 

4. 

4x 

4. 

4. 

/t:  (i  =  x,y,z;  j  =  x,y,z} 


M 


A1A2 

m=x,u 

n=x,u 

M 

M 

xy 

M 

M 

yx 

M 

M 

*-yz 

M 

zx 

M 

^zy 

M 

zz 

=  M-^  (i  =  x,j,z;  j  =  x,j,z)  , 
_sgn(yKj)  y  y  Qi  vj/ 

/  V  /  .  '^mr'-n2  ^  mn 


AT: 


AjA^ 

m=x. 

\  ^ 

,u  n=x,w 

4:, 

■T, 

4. 

4. 

xK, 

•4 

4 

=  yl4j  (i  =  x,j,z;  j  =  x,j,z) 

IeL.e;;'P. 

^1^2  m=x,u  n=x,u 


(271) 


(272) 


(273) 


(274) 


InEqs.  (271)-(274),  'F,„, ,  and  are  the  remaining  integrals  with  respeet  to 


the  parameterized  rough  surfaee  [of.  Eq.  (157)].  The  analysis  above  determines  olosed- 
form  solutions  for  these  integrals  in  sub-seotion  4.3.1  for  smooth-to-moderately  rough 
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surfaces  and  in  sub-section  4.3.2  for  very  rough  surfaces.  With  this  said,  one  determines 
the  other  parameters  in  Eqs.  (271)-(274)  as 


otherwise 

otherwise 


if  A4  =  A4  ,  A4  ,  A4  ,  or  A4 
otherwise 

if  Af-  =  JV  ,  JV  ,  JV  ,  or  Af 
otherwise 


(275) 


1,2 


-I- 


(y-vf +  (z-v)" 


■'xl,!  ’ 


(276) 


A 


x\a 


(y  ■  v)'  (1  +  ^,12 )  +  (z  ■  v)"  (1  +  2 )  -  (1  +  2 )] 

[(y  •  v)  z^,i,2  - (z  •  v)  t?,i,2 ] [(y  •  v)  Z^,l,2  +  (z  ■  v)  1,2 ]  4,: 


^^i,2  =  [(y-v)(z-u)-(y-u)(z-v)] 

(1  +  ^11,2  )  ^U,2  +  (1  +  'ill, 2  )  [(y  •  v)  1,2  -  (z  •  v)  z7^i,2  ]' 
[(y-v)z^,u+(^-’^)^zi,2] 


(277) 


(278) 


A 


xl,2 


-I- 


-I- 


(1  +  'll, 2  )  ^.1,2  [(y  •  v)  t?,i,2  -  (z  ■  v)  z7^i^2  ]' 

{(y  ■  v)  (z  ■  v)  [(y  ■  v)'  (l  +  )  +  (z  ■  v)'  (l  +  )  -  (l  +  2 )]  ^yx,2 

(z  ■  v)"  (1  +  ^,12 )  +  (y  ■  vf  (z  ■  vf  (1  +  ^,12 )  +  (y  ■  v)'  (1  +  V2 )]  '^.1,2}  <2 


,(279) 


'^«i,2=[(y-v)(z-u)-(y-u)(z-v)] 

{[(y  •  V)  (Z  ■  V)  (1  +  /]|i  2  )  Z^,l,2  +  (Z  •  V)"  (1  +  ^11,2  )  -  (1  +  'll,2  )  '^zl,2 

[(y  •  V)  1,2  -  (Z  ■  V)  '^,1,2  ]  -  (y  ■  V)  (1  +  r^i,2  )  '^xl,2}  '^xl,2 


(280) 
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,(281) 


Kl,2=iy-'^y(^  +  ^±l.2)^ylAl,2 

+  (y  ■  v) (z  ■  v) [(y  ■  v)'  (l  +  )  +  (z  ■  v)'  (l  +  ) “ (l  +  ''11,2 )]  4,24,2 

+ [(y  •  v)"  (1 + ) + (y  •  v)'  (z  ■  v)'  (1 + ^,1^2 ) + (z  ■  v)'  (1 + 2 )]  4,24.2 
+  (z  ■  v)'  (1  +  r^i_2 )  4,2  -  2  (y  ■  v)  (z  ■  v)  (1  +  2 )  4,242 

4,2  =  [(y  •v)(z-u)-(y-u)(z-v)]{(z-v)[(l  +  r^i^2)-(y-vf  (1  +  ^2 )]4.2 

+  (z  ■  v)  [(1  +  2 )  4,2  +  (y  •  v)"  (1  +  ?i|i,2 )  4,2  ]  .  (282) 

+  (y  •  v)  [(y  ■  vf  (1  +  ^11 2 )  -  (z  ■  v)'  (1  +  ^11 2 )  -  (1  +  r^i,2 )]  4i,24i,2}  4,2 

4i,2  =  -|(1  -  '2i,2  )  [(y  •  v)  1,2  -  (z  •  v)  4i,2  ]" 

r  n  1  ’ 

+ [(y  ■  v)4  (z  ■  v)'  J  (1  -  ^,12 )  4,2}  [(y  •  y)  4i,2 + (z  •  v)  4,2  ] 

04=[(y-v)(z-u)-(y-u)(z-v)] 

sr  nr  n  ’  (284) 

h,2  -  ''11,2 )  [(y  •  y)  4i,2  -  (z  ■  v)  t?^i,2  ]  [(y  ■  v)  z^^i,2 + (z  ■  v)  1 2  ]  4,2 

0x1,2 = -{(y  •  y)  [(y  •  v)4  (z  ■  v)'  ]  (1  -  ^,12 )  4,2 

+  [(y-vf  (1-^11 2)4,2 +  (z'vf(r^i,2 -V2)  4,2  .  (285) 

- (y  ■  v)  (z  ■  v)  (1  -  2 )  4i,2 ]  [(y  •  v)  4,2  -  (z '  y)^y:,2  ]}  4,2 

04  =  -[(y  •  v)  (z  ■  u)  -  (y  ■  u)  (z  ■  v)]  {(1  -  ^,12 )  4,2  [(y  •  y)^.i,2  -  (z  ■  v)  t?,i,2  ]' 

+  [(y  ■  v)  (z  ■  v)  (?]|1 2  -  /li,2 )  ^yi,2  +  (y  •  v)"  (1  -  ^11,2 )  i?.i,2  ,  (286) 

+  (z'v)'(l-r^i,2)4,2]4,2} 

04  =  {[(y  •  v)4  (z  ■  v)"  ]  (1  -  2 )  [(z  ■  v)  42 + (z  ■  v)  4,2  -  (y  ■  v)  4,24,2  ] 

+  (y  ■  v)  (1  -  r^i_2 ) [(y  •  v)  (z  ■  v)  4,2  -  (y '  v)  (z  ■  v)  4^2  ,  (287) 

+  [(y  ■  v)'  - (z  ■  v)' ]  z7^i,24,2  } 4,2 
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,(288) 


0«1,2  =  -[(y  •  v)  (z  ■  u)  -  (y  ■  u)  (z  ■  v)]  |(l  -  ^,12  )  2  [(y  ■  v)  '  v)  ]" 

+  [(y  ■  v)  (z  ■  v)  (^,12  -  +  (z  ■  v)"  (l  -  ^,1,2 ) 

+  (y-vf  (1-^2)^^!, 2]  4, 2} 

and  ^xu=x-ri,2  ^;,i,2  =  y ' - y  •  v  t?i2  =z-rj_2-z-v.  (289) 

In  Eqs.  (277)-(288),  r^,  2  and  ?jn  2  are  the  perpendieular  and  parallel  Fresnel  refleetion 

eoeffieients  for  two  separate  observation  points,  respeetively.  Refereneing  the  miero- 
seale  seattering  geometry  given  in  Figure  9b,  one  derives  these  eoeffieients  as  [14] 

_  t/eos (^2  1-^0  cos (^/2) 

;7eos(^2)  +  ^ocos(^;2) 


'll!, 2  - 


-//o  eos  {01 2 )  +  0  cos  (0/2 ) 
^oCos(6>'2)  +  /7cos(6>/2) 


-noos(^2)  +  «oA  l-^sin"(^'2) 


‘(^',2)  +  «oA  1-^sin" 


(291) 


where 


and 


cos(^'2)  =  -k;,2-nu  (y-v)^+(z-v)  , 


nn"  (^2)  =  !- cos' (^'2)-!-^^  (y-v)^  +  (z-v) 


(292) 

(293) 

(294) 
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In  writing  Eqs.  (290)  and  (291),  one  must  assume  that  the  medium  below  the  statistieally 
rough  surfaee  S  is  nonmagnetie  with  index  of  refraetion  n  .  Above,  the  medium  is  free 
spaee  with  index  of  refraetion  nQ=\.  One  also  uses  the  law  of  refraetion,  i.e., 

n  sin  (0' )  =  Hq  sin  ) ,  whieh  is  a  direet  result  from  the  geometrieal-opties  approximation 


[22],  Note  that  these  assumptions  are  valid  at  optieal  wavelengths  for  dieleetrie  materials 

[20], 

Specialization  to  conductors 

For  a  very  good  eonduetor,  the  eonduetivity  approaehes  infinity,  a  00  [14]. 

As  sueh,  the  tangential  ineident  eleetrie  field  approximates  to  zero  all  along  the  surfaee  of 
the  eondueting  material,  i.e.,  hxE'  ~  0  ,  and  only  the  eleetrie  equivalent  surfaee  eurrent 
density  J  radiates  when  using  the  PO  approximation.  This  says  that  only  the  dyadie  A/” 
exists,  i.e.,  C~K~  A4  ~  0 ,  and  the  analysis  simplifies  greatly  from  that  of  dieleetries. 
Furthermore,  0^2  ~  0  in  Fqs.  (290)  and  (291)  [14].  Based  on  this  knowledge. 


and 


_  no  cos  (^2) 

-neos^^/j)  no  cos(ei2)-n 
no  eos  ( 2 )  +  cos  ( 2 )  no  eos  ( 2 )  + 

-noos(^j'2)  +  '^o  cos  (0/2 )  -noos(^(2)  +  '^o 
neos^^i'jj  +  ^o  cos  (0/2)  noos(0i'2)  +  ^o 


(295) 


(296) 


21.  Within  a  nonmagnetie  medium  the  permeability  beeomes  that  of  free  spaee,  i.e., 
ju~  jUq,  and  the  index  of  refraetion  n  and  impedanee  r]  simplify,  sueh  that 

n  =  yJeju/leQjUg)  ~  yJs/So  and  rj  =  -sjjije  ~  ^[luJs  (of  Footnote  2,  p.  5)  [14]. 
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Note  that  these  approximations  are  typieally  valid  at  optieal  wavelengths  for  eondueting 
materials.  The  analysis  explores  this  validity  further  in  Appendix  E  below. 
Specialization  to  a  perfect  electrical  conductor 

In  general,  a  PEC  has  infinite  eonduetivity,  and  the  tangential  eleetrie 

field  beeomes  zero  all  along  its  surfaee,  i.e.,  hxE“'  =  0  [14].  Based  on  this  knowledge, 
?]ij  2  =  -1  and  2  =  “1  in  Eqs.  (290)  and  (291).  This  gives  rise  to  a  perfeetly  refleeting 

surfaee,  and  the  analysis  simplifies  even  more  from  that  of  eonduetors.  Eor  a  PEC,  only 
the  dyadie  M  exists  within  the  analysis,  i.e.,  C  =  K  =  AT  =  0 .  Moreover,  the  dyadie 
AT  simplifies  [ef  Eq.  (274)],  sueh  that 


where 


Af: 


A 

=  J/'j  (i  =  x,y,z;  j  =  x,y,z) 

(297) 

_4sgn{X-i}^  ,, 

Q  Q  “ml“n2  mn 

^z\ m=x,u  n=x,u 

|-1  if  Af  =  yl^,  a;,,  a;,,  or  A;^ 

[  1  otherwise 

(298) 

(299) 

[I] 

II 

0 

(300) 

2^i,2  =  (y-v)z^,i,2’ 

(301) 

=  [(y-v)(z-u)-(y-u)(z-v)]42, 

(302) 

(303) 
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and  S:i2=[(y-v)(z-u)-(y-u)(z-v)]z^^i2-  (304) 

Note  that  these  assumptions  are  valid  at  optieal  wavelengths  for  a  perfeetly  eondueting 
material. 
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Appendix  E.  Examining  the  validity  of  the  approximations  used  when 

specializing  to  conductors 

The  purpose  of  this  appendix  is  to  explore  the  validity  of  the  approximations 
made  when  speeializing  to  eonduetors.  For  example,  Figure  24  plots  the  relationships 
found  Eqs.  (295)  and  (296)  for  a  typieal  in-plane  seattering  geometry,  where  r  =  Tj  j , 

(p'  =  90° ,  O'  =  45° ,  and  .064  pm  .  The  eonduetors  studied  in  this  example  are 
titanium  (Ti),  aluminum  (Al),  niekel  (Ni),  and  silver  (Ag),  whieh  have  the  following 
eomplex  indiees  of  refraetion:  rij-  =3.388-73.331,  =1.376-710.21, 

=  1.376 -71 0.21 ,  and  =  0.2342  -  77.214  [143].  Note  that  this  example  maintains 

average  pereentage  errors  mueh  less  than  one  for  all  of  the  eonduetors  studied. 

Figure  25  shows  results  for  a  taetieal  engagement  seenario  with  anisotropie  beam 
parameters,  so  that  A^=A^,  B^^=B^  =  0,  w^=  IMjl  cm  ,  =  0.5w, ,  f  „„  =  0.25w, , 

Aq  =1.064  pm,  =1.376-710.21,  =r^2  =10  km.  O'  =45°,  =\00A^,  and 

<T^  =  lO/lg .  This  setup  eorresponds  with  unpolarized  partially  eoherent  illumination  at 

non-normal  ineidenee  of  a  very  rough  eondueting  surfaee.  As  sueh,  the  approximations 
leading  up  to  Eqs.  (295)  and  (296)  are  typieally  valid  when  eonsidering  the  magnitude  of 
the  seattered  speetral  degree  of  eoherenee  (SDoC)  (Eigure  25a),  the  normalized  seattered 
speetral  density  (SD)  (Eigure  25b),  and  the  seattered  degree  of  polarization  (DoP)  (Eigure 
25e)  with  eonduetors.  The  results  displayed  here  show  exaet  agreement  with  respeet  to 
the  seattered  SDoC  and  SD  but  not  the  seattered  DoP.  Thus,  as  a  rule  of  thumb,  one 
should  use  the  analysis  eontained  in  the  general  problem  to  avoid  errors  in  the  results. 
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Figure  24.  Assessment  of  the  approximations  used  with  the  parallel  and 
perpendicular  Fresnel  reflection  coefficients  for  titanium  (Ti),  aluminum  (Al),  nickel 
(Ni),  and  silver  (Ag).  (a)  depicts  the  amplitude,  whereas  (h)  depicts  the  phase  as  a 
function  of  a  single  polar  angle. 
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0"  [deg] 

(c) 

Figure  25.  Assessment  of  the  approximations  used  when  specializing  to  conductors, 
(a)  shows  the  magnitude  of  the  scattered  spectral  degree  of  coherence  as  a  function 
of  the  difference  between  two  polar  angles,  whereas  (h)  shows  the  scattered 
normalized  spectral  density  and  (c)  shows  the  scattered  degree  of  polarization  as  a 
function  of  a  single  polar  angle.  Note  that  the  analysis  includes  results  for  a  perfect 
electrical  conductor  (PEC)  for  comparison  purposes. 
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